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Abstract 

Light front formalism for composite systems is presented. Derivation of equations for 
bound state and scattering problems are given. Methods of constructing of elastic 
form factors and scattering amplitudes of composite particles are reviewed. Elastic 
form factors in the impulse approximation are calculated. Scattering amplitudes 
for relativistic bound states are constructed. Some model cases for transition ampli- 
tudes are considered. Deep inelastic form factors (structure functions) are expressed 
through light front wave functions. It is shown that taking into account of trans- 
verse motion of partons leads to the violation of Bjorken scaling and structure 
functions become square of transverse momentum dependent. Possible explanation 
of the EMC-effect is given. Problem of light front relativization of wave functions 
of lightest nuclei is considered. Scaling properties of deuteron, ^He and '^He light 
front wave functions are checked in a rather wide energy range. 
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I. Introduction 



Traditional way of describing bound state in quantum field theory is the 
Bethe-Salpeter formalism [1] . However, the dependence of the Bethe-Salpeter 
wave function on the relative time of two particles makes difficult the prob- 
abilistic interpretation of the wave function and gives rise to a number of 
mathematical problems. On the other hand, the relativistic quasipotential ap- 
proach of Logunov-Tavkhelidze [2] based on the two-time formalism for Green 
functions enables one to avoid the difficulty connected with the relative time 
and is much simpler in applications. Similar three dimensional approach for 
bound state and scattering problems in Hamiltonian formalism of quantum 
field theory has been developed in Ref. [3]. 

We give here the derivation of the three-dimensional equations in the light 
front formalism. This form of the three-dimensional approach turned out to 
be effective in the treatment of high energy interactions of composite systems. 
In particular, it gives a number of scale invariant predictions for observable 
quantities which can be checked experimentally. In this approach the relativis- 
tic composite system with the total 4-momentum P is described by means of 
the relativistic wave function ^p{[xi,pi^±]) [1]-[I0], where the "longitudinal 
motion" of constituents is parametrized by means of the scale-invariant vari- 
ables 

_ Pi,0 + Pi,z 

where p^,^ (/i = 0,1,2,3 is the Lorentz index) and are an individual 4- 
momentum of the i-th constituent and the total momentum of the system, 
respectively. Variables Xi are ratios of the light front variables. In terms of these 
variables the wave function of the composite system reflects, in particular, the 
dependence of the internal motion of constituents on the total momentum of 
the system. 

In what follows upper-case letters will denote the characteristics (momenta, 
masses, etc.) of the composite systems and lower-case letters the characteristics 
of constituents. Square brackets in the argument of the wave function $p 
denote the set of the variables Xi and pi^± which satisfy the conditions: 

A A ^ 

J2xi = l; < Xi < 1; J2Pi,^ = ^^- 

i=l i=l 

Recent research has given an insight into the complicated relativistic structure 
of hadrons. Besides, the study of high-energy nucleus interactions has shown 
that idea of nuclei as systems of quasi-independent nonrelativistic nucleons is 
incomplete. Hence it has became necessary to take into account the relativistic 
character of intrinsic motion of nucleons and to consider the quark degrees of 
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freedom in nuclei. 



The above-stated reasons have motivated the subject of research considered 
in this review. 

According to modern concepts nucleons and nuclei are bound states of more 
elementary constituents. Interactions between them or their interactions with 
particles which at the present stage can be considered as elementary are prin- 
cipal sources of information on bound states. 

Investigations with beams of high-energy nuclei carried out in Dubna, Berke- 
ley and other laboratories since the beginning of the seventies have stimulated 
interest in the study of nuclei as relativistic composite systems and the pro- 
cesses of interaction between them (see, e.g., review papers [7|, [0], P^T] 
and the references therein). Presently, there are beams of high-energy nuclei 
with various atomic number in a wide energy range from several GeV/nucleon 
up to several hundred GeV/nucleon at CERN [36j-[38j. Note that this range 
corresponds to relativistic physics and the non-relativistic theory of nuclear 
reactions (see, e.g., the review papers [52], [ID] and references therein) based 
on the Schroedinger formalism of quantum mechanics does not work in this 
case. Relativistic effects have to be taken into account in this energy region 
and high-energy phenomena can be adequately described only using the rela- 
tivistically invariant formalism. It seems more convenient to plot also the data 
in terms of variables which have the relativistic origin. 

In the experiments with beams of high-energy nuclei the situations often 
occur when small internucleon distances must be taken into account. This 
corresponds to large relative momenta of constituents of the nucleus. The 
problem arises to describe adequately nuclei at arbitrary momenta of their 
constituents. That is why it is more convenient to describe the processes in- 
volving relativistic nuclei in terms of wave functions, in which the relativistic 
character of nucleon motion is taken into account, instead of the ordinary 
quantum-mechanical nuclear wave functions, which correspond to the motion 
of nucleons with small internal momenta. 

Note that the idea to use light front formalism in the relativistic nuclear 
physics has been put forward in Ref. [7]. At present there are a number of 
interesting reviews on the light front dynamics and its applications (see, e.g., 
Refs. |ll]-[50]). Some interesting applications of light front variables in the 
high energy hadron-hadron and nucleus-nucleus phenomenology can be found 
in Refs. [51j, i52j and references therein. 

The review is organized as follows: 

Chapter II is devoted to the formulation of the three-dimensional formalism 
for composite systems in terms of light front variables. Equations are derived 
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for boTind states and scattering problems. Cases of spin-0 and spin-1/2 con- 
stituents are considered. It is shown how equations of this approach are related 
with or differ from the Weinberg equation obtained in the framework of the 
old-fashioned perturbation theory in the infinite momentum frame. 

Chapter III deals with the method of constructing relativistic elastic form- 
factors and scattering amplitudes of composite systems in the light front for- 
malism. A general expression is obtained for the matrix element of the current 
of composite system in terms of light front wave functions and the generalized 
vertex operator F^. The explicit form of the vertex operator is found in the 
impulse approximation and the electromagnetic form-factor is calculated for 
a system, consisting of two or arbitrary number of constituents. 

Asymptotic behaviour of the pion form-factor is investigated at large momen- 
tum transfer. 

Problems of the interaction of relativistic composite systems are also discussed 
in this chapter. The scattering amplitude is expressed in general form using 
relativistic wave functions and the transition operator. The constituent inter- 
change mechanism is considered. 

Chapter IV is devoted to the study of deep inelastic from-factors of compos- 
ite systems. Like the case of elastic form-factors, the general expression for 
the deep inelastic tensor W^i, is calculated in terms of the relativistic wave 
functions and the generalized two-photon vertex F^,^. The explicit form of 
this operator is calculated in the lowest order in the electromagnetic inter- 
action and expressions are given for the structure functions Wi and i'W2- It 
is shown that if the transverse motion of quarks is taken into account, the 
Bjorken scaling is violated and the structure functions become the square of 
the momentum transfer dependent. The problem of explaining of the so-called 
EMC-effect is stated. It is shown that one possible way to explain this effect 
is the scattering of electrons on the colourless multiquark configurations in 
nuclei. 

In Chapter V problems of the relativization of nuclear wave functions are posed 
and solved. The case of the simplest nucleus - deuteron is considered in detail. 
A way to the relativization of the well-known nonrelativistic wave functions 
of the deuteron is indicated. The scaling properties of the light front wave 
function of the deuteron are checked by studying spectator-nucleon distribu- 
tions in the deuteron break-up at various incident energies. It is shown that 
the scaling properties of the relativistic wave function hold with good accu- 
racy in a wide energy range. Calculations are performed using the relativistic 
analogue of the well-known Hulten wave function of the deuteron. 

Problems of relativization of wave functions of more complex nuclei are also 
considered. The differential cross section of the break-up of the relativistic nu- 
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cleus on the hydrogen target is calculated. Theoretical calculations arc com- 
pared with experimental data on the "^He '''He. ^He d fragmentation. 
The analysis of data shows that the scaling properties of the relativistic wave 
functions of the ^He and ^He nuclei hold with a good accuracy in the energy 
range considered. 



1 Notations 



Let us introduce the light front coordinates. Instead of usual Lorentz coor- 
dinates Xfj,{t,x,y, z) = {xq,Xi,X2,Xs) we consider new light front coordinates 
), where x^ = ^ {xq ± 0:3). If are the Lorentz coor- 
dinates, then transformation to light front coordinates is performed by the 
matrix C, 



Xii 



— ^ ^ 

10 
1 

V72 -72/ 



Metric tensors are related to each other by the relation 

^1, 0, 0, ^ 
0, -1, 0, 
0, 0, -1, 



9lJ.v — Ca^iQapCal , g, 



\0, 0, 0, -1) 



9nu 



'^O, 0, 0, 1^ 

0, -1, 0, 

0, 0, -1, 

VI, 0, 0, oy 



It is obvious that 



X — — X 

xy = x+y^ + - x±y±. 



2 Generators of the Poincare Group 



Let's consider generators of the Poincare group. In usual notations P^, 
obey the algebra: 
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[M^^, Pp] = i{gr,pPf, - g^pPu), 



Generators of three dimensional rotations and pure Lorentz transformations 
(boosts) are related to Mp,i, as follows: 



Ji = - Sijk Mjk = (M23, M31, M12), 



Oi 



or 



M, 



0, -Xi, -i^2, -i^3^ 
Xl, 0, J3, -J2 

X2, -Ja, 0, Ji 

l^Ks, J2, -Ji, J 



Performing the transformation by means of C-matrix, we obtain for the light 
front generators: 

P^={P^,Pr,P2,P-), P±^^{Po±P3), 

( 0, -Sr, -S2, Ks^ 
Si, 0, J3, Bi 
0, B2 



where 



M, 



Bi^ -^{Mio + Msi), 



I.e. 



^1 = ^(^1 + ^2), S2 = -^(X2- Ji), 
^1 = ^(^1-^2), S2^^{K2 + Jl). 
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In the usual covariant theory initial state of the system is given on the hyper- 
plane t = const. In quantum mechanics and quantum field theory canonical 
commutation relations between canonically conjugated quantities are given 
on this hyperplane. The plane t = const remains invariant under three- 
transformations generated by momentum operators Pi and under three di- 
mensional rotations generated by angular momentum operators Mij = —Mji, 
e.i. the stability group of the hyperplane t = const is 0(3) x T(3) (or SU(2) x 
T(3)). 

Corresponding operators and Pi are called kinematic operators (according 
to Dirac). If the system of particles is considered, eigenvalues of kinemati- 
cal operators are equal to the sum of eigenvalues of kinematical operators of 
constituents. 

Evolution of the system proceeds in the direction orthogonal to the initial 
hyperplane, or in the direction of the time. The operator, which generates this 
evolution is Hamiltonian H = Pq, H is dynamical operator. For the system of 
particles it is not equal to the sum of Hamiltonians of individual particles, but 
depends on the interaction between them. Besides Pq Lorentz boost operators 
Ki are also dynamical operators. So we have 6 kinematical operators {Pi, Mij) 
and 4 dynamical operators {Po,Ki). 

Another possibility for postulating commutation relations was pointed out 
by Dirac in 1949 [53]. In particular, he proposed to use light front hyperplane 
x+ = const for postulating commutation relations. The hyperplane a;+ = const 
is not changed under translations in the orthogonal directions x±, which are 
generated by P± and in x_ direction, which is generated by P+. Besides it is 
invariant under two-dimensional rotations generated by M12 = J3 and boosts 
generated by Bi {i = 1,2). Corresponding stability group is E{2) x T(3). 
Evolution of the system proceeds in the x+-direction and the role of shift 
generators in this direction is played by P~ = P+. Thus we have 6 kinematic 
operators {P~^, P±, J3, Bi) and 4 dynamical operators {P~ , Si, K3). 

From commutation relations between generators: 

[p-,P''] = [p-,P+] = [J,,p-]=0, 
[B„p-] = -iP-, 

it follows that an arbitrary 4- vector {^'^ ■, ^~) is transformed according 
to the rule: 

It is obvious that B12 generate Gahlei boosts in the directions x and y, P+ 
plays the role of mass, plays the role of momentum, J3 plays the role of 
angular momentum and P~ plays the role of the Hamiltonian. 



8 



Consider now commutation relations between dynamical operators: 



[s„s,] = o, [s„p-] = o, 

[J3, Si] = iSijkSk, 

[5„p,] = -^%p+. 

From these commutation relations it follows that transformation rule for ar- 
bitrary 4- vector A'^ is the following: 

Comparison shows that the system of operators (P_, P\ J3, Si) allows the same 
interpretation as {P^, P*, J3, Bi), but in relation to the hyperplane x~ — const, 
i.e. if the initial conditions are given on the plane x~ = const. 

Let us consider now commutation relations: 

[Si, Bj] = -iSijsJs + idijK^, 
[K,,P^] = [K,,h]^Q, 
[K:,,P^]=^iP^. 

Prom these relations it follows that the transformation rule for arbitrary 4- 
vector A'^ looks as follows: 

Besides, from commutation relations: 

[K^,Bi\ = -iBi, 
[K^,Si]^iSi 

it follows that 

e^'^^^BiC-^^'' = e'"Bi, 

Thus, in the light front system the boost operator along the third axis acts as 
a scale transformation. This means that the light front is well suited for the 
treatment of high energy problems. 
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3 Light Front Quantization 



For simplicity let's consider first the free scalar field with mass m. In covariant 
formulation the commutator of fields looks as follows: 

Mx),v{y)] = - y), 

where [M] : 

A(x) = —e{xo)5{x^) ^e{x,)e{x^)Um^). (1.1) 

Using formulae: 

e{xo)e{x^)\^^=Q = 0, e{xo)6{x'^)\^^=o = '^£{x-)6{x^), 
one obtains from (11.11) 

[v{x),v{y)]x+=y+ = -^eix- -y-)S^^\x^-y^). (1.2) 



As it can be seen from here the commutator does not vanish at x_ ±oo. 
Thus, the problem of boundary conditions at infinity arises. It was shown 
[55] that the limit x_ ±oo corresponds to the contribution of zero modes 
P+ = 0. It was shown in Refs. [SB], [5Z] that the commutator on the light 
front depends on the interaction. Various methods of excluding zero modes 
based on the periodic boundary conditions have been proposed. Some authors 
assume that it is necessary to work with two Hamiltonians P_ and P+, or, if 
one uses only one Hamiltonian, say P_, it is necessary to require fulfilment of 
periodic boundary conditions with respect to the variable It is interesting 
to note that if zero modes P+ = are excluded, then the vacuum, as a state 
with P+ = 0, is at the same time the lowest state of the Hamiltonian P_. 

The scattering matrix in the light front quantization is defined as a general- 
ization of Dyson formula: 



S = exp 



— i y d^xHiix) 



:i.3) 



where Hj{x) is the the interaction Hamiltonian in the interaction represen- 
tation, T+ is the ordering operator with respect to the variable x+. If the 
interaction Lagrangian does not contain derivatives, it can be shown that 
Hi{x) = —Ci{x) as in the usual case. For the equivalence with usual case it 
is necessary to show equivalence of T+-product to the usual T-product. 
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Consider 



T{ip{x)(p{y)) = e{xo - yo) Mx), (p{y)] + ^{y)^{x). 

These two expressions may differ only in the cases: when xq > yo, we have 
x+ < I/+ or when xq < yo, we have x+ > y+- It is easy to see that in both cases 
(x — y)"^ < 0, i.e. interval is space-like and commutator is identically zero. 

Thus, both orderings coincide and the light front 5'-matrix coincides with the 
covariant S'-matrix. In the spin case additional noncovariant terms arise, but 
they cancel with corresponding noncovariant terms in propagators [58] |62j. 



II. Light Front Formalism for Bound States 

1 Equation for the Two-Body Bound State Wave Function 

In this section we develop the three-dimensional formalism for composite sys- 
tems in light front variables. These variables were introduced by Dirac [53] 
with the aim to construct the quantum theory with canonical communication 
relations on the light front hyperplane (instead of traditionally used t = hy- 
perplane). We note that the real progress in this direction has been achieved 
much later [5S] [53]. 

Let us consider the Bethe-Salpeter amplitude (wave function) 

XP,a = (0 I T{ip,{x,)ip2ix2)) I P, «) = e-*^^xp,a(^)- (2.1) 

Here | P, a) is the state vector of two particles with total 4-momentum P 
and quantum numbers a, X = {xi + X2)/2 is the centre-of-mass coordinate, 
P — Pi + P2. Relative coordinate and momentum are defined as 

x = xi-X2, p = — - — , (2.2) 
and the light front variables can be introduced 

^^^xo±xs^ P±=Po±P3, P± = Po±A. (2.3) 

Let us introduce the Fourier transform XP,a{p) = XP,a{p-,P+,P±) of the 



11 



Bethe-Salpeter amplitude 



dp^ dp_ dp± e 



-i{p+x-+p-x+-p±xj_) 



XpAp) (2-4) 



and define the function 

/oo 
dp-XpAP-^P+^P±)- 
-oo 



(2.5) 



It can be sliown tliat tlie function pAp+^pA depends on tlie values of the 
Bethe-Salpeter amplitude on the light front hyperplane 

Xo+ X3 = 0. 

In fact, using definition (I2.5P and Fourier transformation (12 ■4p . we get: 



^P,a(P+,P±y 



(27r)3 



dx+ dx- dx± 6{x+)e 



Xp,a(x+,x_,a;_L). 

(2.6) 



Let us consider now the two-particle Green function 



G{xi,X2;x[,X2) = G{X - X';x,x) = (o \ T{ipi{xi)ip2{x2)vi{x[)ip2{x'2)) I 0^ 
= J dPdpdp'e-'''^''-'''^-'^"''''^G{P;p,p'). (2.7) 

Here the total and relative 4-momenta and 4-coordinates in the initial and 
final states are introduced as follows: 

„ , Pi -P2 y X1+X2 , . 

P = Pl+P2, P= , X= , X = Xi-X2, (2.8) 

P =Pi+P2^ P= — 7, — , X= — - — , X =x^-X2. (2.9) 



Fourier transform of the "two-time" Green function can be defined as 

G{P-p+,p^-p'jA= / dp^dpiGiP;p,p'). (2.10) 



For free particles we have 



p,p 



-5^^\p~p') 



i+p 



mi + is 



—p) —m2 + le 



[2.11] 
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Performing the integration according to the definition fl2.10p . we obtain [1]: 



P^x(l — x] 



G^'\P;p+,P±)6{p+-p'^)6^'\p±-p'^). 



(2.12) 



In this expression the variable x is introduced in the following way 

1 , P+ 
^=2 + /^- 

It is obvious that when the variable x varies in the limits 

< X < 1, 

the variable p+ varies in the interval [— P+/2, -P+/2]. 
Inverse operator can be defined by the relation 

I-pJ^'^ I dplG-\P;p^,pr,pl,f[)G{P;pl,p'[-,p'_,,p'^ 



(2.13) 



(2.14) 



5{p+-p'+)5^'\p±-p'a.)- 



(2.15) 



If we introduce the interaction kernel V (quasipotential): 

G-\P;p+,pr,p'^,p'^) = &'^'\P;p+,p±)6{p^~p'^)6^'\p±-p'^) 

-^^V{P;p+,pr,p'+,p'i.) (2.16) 



after simple transformations the equation for the wave function 
$P,„(x,pl) = P+x(l - x)^p_„(p+,p_l) 

takes the form 



(2.17) 



p2 



(pl + (1/2 - x)Px)' + ml (pi + (1/2 - x)P^Y + ml 



X 

r^^TT^ Jdp±V{P;p+,p±;p'^,pl)<l>pA^',pl). (2.18) 



1 — X 



$P,a(a;,RL) 



The equation obtained gives the wave function of a bound state in an arbitrary 
Lorentz reference frame. Comparing it with the equation in the frame where 
P_j_ = we get the transformation property for the wave function from the 
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arbitrary frame to the frame in which the total transverse momentum of two- 
particle bound state is equal to zero: 

^p{x,p^) = $p^=o(^'P± + (1/2 - x)P^). (2.19) 



2 Equation for the Scattering Amplitude and Relation to the Equa- 
tion in the Infinite Momentum Frame 

Let us derive now the equation for the two-body scattering amplitude. The 
definition of the scattering amplitude T(P; p, p') in the 4-dimensional covariant 
Bethe-Salpeter formalism looks as follows: 

+ j dydy'G^^\p-py)T{p-p''y')G^^\p-p'\p') 

= G^") (P; {p - p') + G^") (P; p)T{P- p, pOG^") (P, p'). (2.20) 

We define the quantity T {P; p^, p±; p'^, p'j_) by the similar expression 

G(P;p+,px;p'+,p1) = &'\P;P+,P±)S{p+-p'+)6'~'\p^-p'^) 

+ G'-'\P-p+,P^)f{P-p+,pr,p'^,p'^)&'\P;p'^,p'^). (2.21) 

Integrating fl2.20p according to fl2.10p . we get: 

G = G(°^ + G(o)tS(o). (2.22) 
Comparing formulas (12.221) and (I2.2ip . we obtain: 

T = G(°) ■ GmX^^^") ■ G(°)"' . (2.23) 
It can be shown that on mass-shell the following equality holds: 

f = T. 

Let us derive now the equation for the amplitude T. Using the definition 
fl2.16p , we get the equation for the Fourier transform of the "two-time" Green 
function 

G = G(o) +G(o)1/G. (2.24) 

In fl2.24p the multiplication is understood as a three-dimensional integration 
over the corresponding variables x and p±. Comparing (I2.24p and (12.211) one 
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can see that 

f&^^ = VG (2.25) 
from which the equation for scattering amphtude T follows: 

f = V + VG^^^f. (2.26) 

In the frame, where total transverse momentum is zero P± = 0, the equation 
f l2.26p looks as follows: 



T{P;x,p±;x',p'j^) = V{P;x,p±;x',p'^ 
x"il-x'') 



dx" r „ VjP; x,pr, x",p'[)T{P; x",p'[; x\p'^) 

x" 1—x" 



In a number of papers (see, e.g., [6l]-[67]) the composite systems have been 
described on the basis of the called old-fashioned three-dimensional perturba- 
tion theory in the infinite momentum frame, which has been used by Weinberg 
[68] in the relativistic quantum field theory. The equation fl2.27p is the close 
analogue of the equation derived in [68]. In the lowest order or perturbation 
theory equation (12.27^ reproduces the equation from [68] and at the same 
time contains the regular method of constructing the interaction kernel in the 
higher orders of perturbation theory. We will not discuss this point here, but 
recall that as in the canonical three-dimensional approach |2] there exist two 
methods of constructing of the interaction kernel (by means of the "two-time" 
Green function and by means of the scattering amplitude on the mass-shell). 

We note, however, that there exist one substantial difference between the 
equation derived here and equation of [M]- In the light front approach the 
equation is written in an arbitrary Lorentz frame and "longitudinal motion" of 
constituents is parameterized in terms of scale invariant and Lorentz invariant 
variable x = (-P/2 + p)_|_/P+. The Weinberg equation is written in the infinite 
momentum frame and "longitudinal motion" is parameterized in terms of the 
variable x = {P/2 + p)-i/ P^, which is not Lorentz invariant. 



3 The Case of Two Spin-1/2 Particles 



Let us consider the problem of two spin-1/2 particles interacting with scalar 
or vector gluon fields. 

Let us start with definition of Green function in the interaction representation 
in terms of a;+-ordered product: 

G+{xi,X2] X^.Xa) 



(0 


T+{'lpl{Xi)lp2{x2)llj2{.X^)'^^{Xi)S) 


0) 


(0 


s 


0) 
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where S*- matrix is defined by the relation (11.30 . It was noticed that in models 
like gipTilap, gip'^^ipB^, etc. noncovariant expressions in vertices and propaga- 
tors cancel in S'-matrix. That is why we can use the covariant S'-operator. 

Using usual consideration |1] one can derive the light front Bethe-Salpeter 
equation: 

G+{xi,X2]X^,Xi) = S^p\xi -X3)SP{X2 -X4) 

-Jdx^ dxe dxj dxg S p\xi-x;i)S p\x2-xe)K{x5, xe, X7, xs)G+{x7, Xs;x3, X4), 

where the kernel K contains the sum of all two particle irreducible diagrams 
in the light front formalism. 

Sf{x) = Sf{x) — ^'yJ^e{x^)5{x-^-)6^'^\x±). 

In the momentum representation the free fermion propagator is of the form: 

-n r \ _ P + m l7+_ f + m 

^f[P) — -5 ^r~~~- — o — ~ ~2 2 , • ' 

pz _ _|_ ^ p+ — + te 

where in the nominator the momentum p is on the mass shell: 

^1 1 ^ ^ p i + 

2 2 p+ 

In the momentum representation the Bethe-Salpeter equation looks as follows: 

G+(P;p, q) = S?{f,,P + p)Sf{^i2P - q)S^'\p - q) 

-S?ifi,P + p)SPifi2P-p) J d\'KiP;p,q')G4P;q',q). 

Here the usual notations are used: 

rrii 



P = Pi+P2] P = /i2Pi - /iiP2; /ii + yU2 = 1; /i 



* mi + 1712 



Solving this equation by iterations we get the resolvent representation: 

G^{P-p, q) = 'SfiliiP + p)Sf{lJi2P - q)S^'\p - q) 

- S^ifiiP + p)SP{fi2P - p)T+{P;p, q)S^P{fiiP + q)SP{fi2P - q). 

Here the resolvent T+ coincides with the covariant T-matrix. Thus, the only 
difference consists in the external propagators. 
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The contribution of the two-body bound state can be evaluated in the usual 
way. Let Xi+,X2+ > X3+,X4+. Then: 

P,a 

where bound state wave function is defined as follows: 

Xp,a(xi,X2) = (0 I T+(V^(xi)^(x2)) I P,a). 
For the contribution of bound state Mb we have: 

xe( x,-x' 



1 








2 F+ 





Here usual notation of centre of mass coordinates (X, X') and relative coor- 
dinates (x, x') are introduced. This expression can be rewritten as: 



J d''P9{P+)e 



+.rp(x-x') XpA^JXpA^ j ^■(p-_p-)(|^^|+|^vi) 



{2tt)^J ' ' P+{P- - Ps +ie) 



Pl + Ml 



Thus in the momentum representation we have: 



Equation for bound state wave function looks as follows [6]: 

XpAp) = -S?{fiiP + p)SP{fi2P-p) J d^qK{P;p,q)xpAQ)- 

Let us proceed now to the three dimensional (quasipotential) light front for- 
mulation of the bound state problem. 

Bethe-Salpeter amplitude (wave function) is defined as: 

XpA^u^2) = (0 I T+{i^{xi)Ax2)) \P,a) = e'^'^xpAA- 



Let us derive the equation for the three-dimensional wave function Xp.a 
0, a;), x = {x-,x±). 
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In the momentum representation we have: 

/oo 
Xp{p)dp- 
-oo 

It is necessary to investigate properties of the "two-time" Green function: 

~ /"OO 

G+{P;p,q)= / G+{P;p,q)dp_dq_. 

J — CO 

Let us consider first the free Green function: 

GoiP;p,q) = (Pi + mi){% + m2)5^^\p - q)Io, 

where 

/o= r dp^dq^ A«(/iiP + p)A(')(//iP-p)5(p_-?_) 

J—oo 

1 r°° dp_ 



L 



XP+(1 - X)P+ 7-00 (p_ -i^,P_-p- + f - IJ^P_ +p- + ^) ' 



x = , 1 — a; = , < a; < 1. 



Special property of the spin case is that numerators of free propagators do 
not participate in the integration, because they do not depend on p_ and g_. 
After integration we obtain: 

Go{P;p,q) = {Pi + mi){p2 + m2) - - 



x{l - x)P+ (p2 _ _ ^ ^g-j ■ 



The reference frame with Pi = is used in this integration. 

Since (7+ (<5_|_) is defined on the subspace with positive frequencies, it is nec- 
essary to project onto this subspace. Define: 

«^HPl)«2Hp2)T«f(9l)«2^(92) = r^^^-^i^^ 

and obtain the matrix equation 

g+ = goxf X go, 

where 

-2Ki 2/;/i2///2^(.r)^(l - .r) 



g = go{P;x,pj_) 



x{l - x)P+ (p2 _ pj+'^l _ pj+ml . X 

V X 1—x ' 
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q, > 1 . P2 

Fig. 1. 

Inverting the operator we get: 

9+ = - V. 

Here the quasipotential V is represented by the series, which can be written 
in the following formal form: 

v = f^{i + ~g,fy' = (1 + fgoY' X r. 

The equation for the bound state wave function is of the form [6]: 
V X 1 — X 



2mi 2m2 



2(27r)3a;(l - x) 



/ dx' / dp'j^ Vp {x,p±;x', p'^ ) $ (x' , p'^ 

1 ~ X) JQ J 



The wave function $p is related to the initial wave function if) by the relation 

^p{.x,pi_) =Ui{Pi)u2iP2)i'pix,P±)- 

This equation coincides with corresponding equation for scalar particles [1]. 

Let us construct now one meson exchange quasipotential in the theory gtpipip, 
which corresponds to the diagram of Fig. 1 and looks as follows: 

K 



(27r)4 {p - g)2 _ ^2 + • 



According to the quasipotential constructing rule it is necessary to consider 
the integral: 

h= r dp.dq. A^o\f^iP + P)^o\f^2P-p)K{P;p,q)S{p--q-) 

J — CO 

xA(^)(/xiP + g)A(')(/i2P-g), (2.28) 

where 

-ig^ Ui{Pi)u2{p2)ui{qi)u2{q2) 



K{P;p,q) 
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Integrating the expression fl2.28p we obtain for the quasipotentiah 

Vp = g^Ui{p^)u^{q;)U2{p^)u2{qM^)G{'^ - x)e{x')e{l - x')Vp, 

where 



9{x — x) 



p2_ 



\x — x'|V(x, p_L,x',p'j_) 
1 — X 1 — x' 



+ ie 



+6(x' — x) 



X — X 
I \2 I ,,2 



p{ + mj p!l + ml {p± - p'^Y + jj, 



x' 



X' — X 



+ ie 



. (2.29) 



Note that hght front formulation of the quasipotential approach for spin par- 
ticles is free from difficulties, which are characteristic to the usual formulation 
with equal times [69] . 



4 Equation for the Many-Body Bound State Wave Function 



Formalism developed in the previous sections can be generalized to the case of 
N relativistic interacting particles. The way of this generalization can be seen, 
if instead of the variable x, defined by the relative momentum, two variables 



x^^^ and x^'^\ defined by the individual momentum of particles are used: 



(0 
-t -I- 



(2.30) 



The variables x*-*-* vary in the interval < x*-*^ < 1. 



Define the Fourier transform of many-body Bethe-Salpeter amplitude (wave 
function) : 



XP,.([x«]) = (O I T(v.i(x«)v.2(a:f ) ■ ■ ■ ^v(xW)) | P, a 
by the following equation 



N 



n d^x^'^ exp 



N 

J2p 



I 

L i=l 



XP,a([x«]), (2.31) 



where 



[p»]=p«,...,p(^); [x(:)]=x«,...,x(f). 
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Here we have ascribed the Lorentz index to the 4-coordinates x*^*) in or- 
der to distinguish them from the scale-invariant variables x^^\ which will be 
introduced later. 



If we introduce the light front variables 



V'') I 

X — — 'T' 
n — ^ 



(2.32) 



N 



and integrate (12.311) over H dp^- ■ We obtain 

i=l 



N 



i=l 



N 



i=l 



X_ p j_ X j_ ^ 



XpA{4^])- (2.33) 



The function \E'pQ,([p^'',p|,^^]) is related to the Bethe-Salpeter amplitude in 
the following way 



N 



(i); 



N 



.(0 



(2.34) 



i=l 



i=l 



Let us introduce now the Fourier transform of the "two-time" Green function 

GiP;[p?,p?];[p?\p?']) 

Y[dp^iUp^i^'s(p.-Y,p'-!^]6ip--j2p^:^' ]g{p; [p«]; [p«']). (2.35) 



i=l 



1=1 



i=l 



The function G{P; [p''*-']; [p*-*^ ]) is defined by the Fourier transformation 
G([a:«]; [xf]) = (0 I T(^i(xW) ■ ■ ■ ^^(xf)), ^+(x«') ■ ■ ■ ^^(xf)')) | 



N 



= {2n)-^^ I n t^V'^f^V'^' exp 

i=l 

xG(P;[p«];[p«']). 
For the case of free particles we have 



N 



(2.36) 



— i ^(p'"'x'"' — p X 

i=l 



■N 



N 



G(°)(P;[p«];[p»']) 



^9 9 

n (p(*) - m» + ie) 



(2.37) 
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Integrating both sides of (12.371) according to the definition fl2.35p and omitting 
the 5- function corresponding to the total 4- momentum conservation, we get [5j: 



1=1 



i=l 



N 



p^-' n 

i=l 



P2 _ g ip^-x(^)p^r+mj 
1=1 



&'\P; \p?,P?])I[S{p? -pf)6^'\pf-p?'). (2.38) 



N 



1=1 



The variables x'-*-* are defined in the following way 



X 



1 = 1,2,. ..,N. 



(2.39) 



Thus, the function G''^^{P; [p^+,p^±]) is defined under the conditions 



N 



N 



5]a;(*) = l; < < 1; T.P? = P- 



(2.40) 



i=l 



i=l 



Let us introduce now the inverse operator G ^ by means of the relation 



p. N 



i=l 



N 



Udp?" / Y[dp^'y'G-\P- [p?,p?]; [pf',pf]) 



i=l 



N 



G{P- [pf,pf"]; [pf ,pf ]) = n^(P? -pPV^'Hp? -Pf ) (2.41) 



i=l 



and define the interaction kernel V: 

G-\P; [p?,P?]; [P?',P?']) = -\P- [pf ,pf ]) 



y(P;[p«,p?];bf (2.42) 



(2i)^(27ri)^-i 



The equation for the wave function 

$p,a(k«,P?]) = p^'U^^'^^pAIp^Ip?]) 

looks as follows 



N 



i=l 



(2.43) 
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i=l 



^pA[^^'\p?]) 

I n 1^* 1 - E-'"' / n*i''* - Erf'' 

X V{P- [p?,p?]; [pf,pf ])$P,.([x«',p«']). 



(2.44) 



Corresponding equation for the case of spin- 1/2 constituents looks as follows 



p2 £.(p?-xWPx)^ + mf 



1=1 



(i) 



N 



(0 



J i=l 



1 = 1 \ 1 = 1 / 4 = 1 \ 1 = 1 / 



N 



N 



j=l i=l 



The formalism developed can be used for the treatment of a wide class of 
elementary particle and nuclear physics problems (see, e.g., [72]-[7i]). 

Other forms of the light front description of particle dynamics can be found 
in Refs. [TS]-[7S]. 



III. Relativistic Elastic Form Factors and Scattering Am- 
plitudes for Composite systems 

1 Formulation of the Method 



The reaction of a system on a weak external perturbation corresponding to 
the local field A{x) is described in quantum field theory by the expression 



5S 
5A{k) 



{27T)'6^^\P-Q-k){P,a I J(0) I Q,P). (3.1) 



A=0 



Here J{x) is the local current of the system 

. . . SS 
J[x) 



5A{x) 



(3.2) 



|P, a) and are the state vectors of composite particles with momenta 

P and Q and the sets of additional quantum numbers a and respectively. 
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normalized in a relativistically invariant manner 

{P,a I Q,/3) = 2Po(27r)35(3)(p_Q)5^^. 



(3.3) 



Below we suggest a method of constructing relativistically covariant form fac- 
tors of composite systems in terms of "equal time" three-dimensional wave 
function ^p^a{[x''^\Pi^])- 

Let us consider first the case of two-particle system and introduce the quantity 
R defined by the vacuum expectation value of the chronologically ordered 
product of Heisenberg field operators of the scalar particles (pi{xi) and a some 
local current J{x): 

R{xi,X2;yi,y2) = (0 | T{(pi{xi)(p2ix2)JiO)(pt{yi)(ptiy2)) \ 0) 

= (27r) / d'^pi d'^p2 d'^qi d'^q2 exp — i'^iPjXj — qjy^ 

I j=i 

X R{pi,p2]qi,q2). (3.4) 



Introducing, as above, the relative 4-coordinates and 4-momenta 

y X1 + X2 ^ yi + y2 

X = ^ , X = Xi-X2, Y = , y = y^-y2] 

„ , P1-P2 ^ . qi-q2 ^ ^ 

P = vi+V2, p = — - — , Q = qi + q2, q = — - — , 

we rewrite expression fl3.4p in the form 

R{X,x;Y,y) 

= {2n)-^^ fd^Pd'^pd^QdS ew \ - i{Px-QY+px-qy)]R{P,p;Q,q). (3.6) 



As is known [79j, the R quantity can be presented in the form 

R = GTG, 

or in the detailed form: 

R{X, X] Y,y)= J d^X' d^x' d^Y' d^' G{X - X'; x, x') 
xT{X',x';Y'y')G{Y'-Y;y,y'). 



(3.7) 



(3.8) 



In the momentum space we get 

R{P,p;Q,q) = J d^p'd'q'G{P-p,p')T{P,p';Q,q')G{Q-q',q). (3.9) 
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Fig. 2. 

Here G is the two-particle Green function of scalar fields (pi{xi) and the vertex 
function F is the sum of all two-particle irreducible diagrams for 5-point Green 
function (see Fig. [2]). 

Passing to the "two-time" description in terms of light front variables, we 
define the quantity 

~ r+oo 

R{P;p+,p±]Q]q+,q±) = dp.dq. R{P,p]Q,q). (3.10) 

J —oo 



The quantity R can be presented in the form 

R = GfG (3.11) 

or in the detailed form 

/■-P+/2 r rQ+/2 r 

RiP;p+,p±;Q;q+,q±) = dp'+ dp'i. / dp'l dpi 

X G{P]p+,p±]p'^,f^)f{P]p^,p^]pl,p[)G{Q-,pl,p[-,q+,q±). (3.12) 

Here F is the vertex integral operator. Let us show that the quantity F defines 
the form factor of composite system. Starting from the spectral properties ^SU] 
of the 5-point Green function fl3.10p . it is possible to show that the quantity 
R has the pole singularities near the points corresponding to the masses Ma 
and Mf3 of composite systems: 

R{P;p+,p±;Q;q+,q±) 

r... ,412 ^P,a(p+,P±)(^,c^ I ^(0) I Q , (3)% ^p{q+, qA_) 
^^'^ ^^^ {P' - Ml){Q^ - Mj) ■ ^^-^^^ 

On the other hand, taking into account the pole singularities of the two- 
particle "two-time" Green function 

G(P;p,.p,;p;.pl).^^^<2±fM>:^. (3.14) 
we find from (13.121) 

R{P]p+,p±]Q,q+,q±) 
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,412 ^P,-(P+'P±)^Q,,3(?+,9l) r^+Z^ r ^^ rQ+/2 



x'i't„(p'+,p'±)r„e(/';p'+.p'x; 0; 9V, fD^Q.sl^i, fl), (3.15) 

where 



ra/3(i^;p+,P±;Q;g+,g±) = r(P;p+,_pl;Q;g+,gl^ 



(3.16) 



Comparing equations (13.131) and (I3.15P we get the following expression for the 
matrix element of the local current J: 

{P,a I J(0) I g,/3) = / / dp^ / / du 

J~P+/2 J J-Q+12 J 

X ^p,a(p+5P±)fa/3(P;p+,p±;<5;g+,'?l)^Q,/3(g+,gl)- (3.17) 

Equations (I3.15P and (I3.17P give an exact expression for the vertex operator 
of the composite system in terms of 4- and 5-point Green functions G and F 

^af3iP;p+,p±;Q;q+,q±) = hm / dp'^ dp'^ / dq'^ dq'^ 

Q'^'^Mp2 

xG-\P;p+,pr,p'^,p'^)[GYV]iP;p'^,p'^;Q;q'^,q'±) 

xG-\Q;q'_^,q'^;q+,q^) (3.18) 

Using the perturbation theory methods for these functions one can construct 
the coupling constant expansion for the vertex function of composite system. 

In the case of spin-1/2 constituents the matrix element of the local cyrrent 
operator looks as follows [71j : 

(27r)2(Q,/3| J^(0)|P,a) 
= / dq_dj^^^,{q) '-1±^ f ,(Q, g; P, pf-l±^ ^p,„(p), 

where, for instance, dq = dq+ dq±, q = {q+,q±). 



2 Elastic Form Factor in the Impulse Approximation 



Let us consider the so-called impulse approximation for the vertex operator T 
which corresponds to the limit of "weakly bound" (noninteracting) particles 
(see Fig.E]). 
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Fig. 3. 

For the vertex operator, corresponding to the conserved vector current, we 
find 

f, = ffj+fg, (3.19) 
ffj = [G(°)]-^[G(o)f£)G{o)] [G^W (3.20) 



where 
Here 



= {2n)\{p, + q.),5^'\p, - q,)[Gf\p,)]-\^. (3.21) 



G(°)(p,) = Gf\p,)Gf\p,) = Ylip' - ^)-\ (3-22) 
G^'\q,) = G^?{q,)GP{q,) = \{{q^ - m^)-' (3.23) 

i=l 

are the two-particle Green function for free particles with masses rrii and 
charges Cj. 

Then for the invariant form vector of the composite system defined by the 
relation 

{P,a\J,iO)\Q,P) = iP + Q),FiA^), A = P-Q, (3.24) 
in the reference frame, in which 

P+ = Q+ , (P - QY = A' = -Al = -{P^ - Q^f, (3.25) 

we have [Ij: 



-A 



2, ei(27r)^ dx 



2 Jo x(l — x) 
X J dp±^p^^Q{x,p± + {I- x)A^)<^p^^q{x,p±) 

+ similar term with 62- (3.26) 



Note that construction of relativistic form factors of composite systems in 
other versions of relativistic description of bound states is considered in Refs. 

ca, m-[ 
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3 Relativistic Form Factor for the Many-Body System 



Let us construct now a form factor for the relativistic many-body system in 
terms of hght front many-body wave functions ^p{[x^^\p^l^]). 

Consider, as in the case of two constituents, the quantity R, which is defined 
by the vacuum expectation value of the chronologically ordered product of the 
Heisenberg field operators (pi{x^^^) and a local current J(x) 



N 

(27r)-^^ / n^^V'^ d^q^'^ exp 

i=l 



N 



=1 



xi?([p«],][g«]). (3.27) 



The quantity R can be presented in the form 

R = GTG. (3.28) 

Multiplication in (13.281) has to be understood as an integration over the 4- 
coordinates of particles. G is the many-body Green function of the fields 
ipi{x^^^) and the vertex function T is defined by the sum of the irreducible 
diagrams of the (2A^ + l)-point function (13.271) . 

Proceeding now to the light front description, we introduce the quantity 
]j [q+Kqa^]) by the relation 

fOO ^ ( ^ \ ( ^ \ 

= l^Iidp^- dq^^sl^P^ - Y.P^:^y[Q- -T.<l^^jR{[p% [g«]) (3.29) 
and write it in the form 

R = GTG. (3.30) 
Multiplication in Eq. (I3.30p has to be understood in the operator sense: 



ab = n dir^Q- - E «?') / n i^fH'^^ - 1:? ?') 

^ 1=1 \ i=l / 1=1 \ 1=1 / 

X A([p?,p?]; [gf ,gf])E([gf ,gf ]; [g«,g«]). (3.31) 



From the spectral properties of the function G [80] it follows that R possesses 
the double pole singularities 
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^ ^'^^ ^ (P^-M^)(Q2_M|) ^^-^^^ 

in the vicinity of the points, where A^-particle system forms the bound states 
with masses Ma and Mp and a set of other quantum number a and /?, respec- 
tively. 

On the other hand, knowing the pole singularities of the Green function 

GibV ,P±]; b+%?V ) -2(27r) — -j- , (3.33) 

p2^M2 [P^ - Mi) 

one can reduce the Eq. (I3.32p to the form 

n*?'* f Px - f:??') / n f - hf] 

•'^ i=i \ 1=1 / •' i=i \ 1=1 / 

•'^ i=i \ i=i I ■' i=i \ i=i ) 

X ^J,J[j>f ,pf ])r.,(bf ,pf ]; [gf ,gf ])*Q.,([gf ,gf ]). (3.34) 

Here 

F + / - \ / / / ■\l \ I / - xl , _ , / — » 



X 



j=l \ i=l / i=l \ i=\ 

^-i^rJ^) ^Wi. rJ^)' 



xG-^([p^^pT];[p^^pT]) 
X [GrG]([pf ,pf ]; [gf ,gf ])G-^([gf ,gf ]; [g?,g?]). (3.35) 



Comparing (13.321) with (13.341) . we get the following expression for the matrix 
element of the bound state current: 

(P,a| J(0)|Q,/3) 

N ( N \ N ( N 



X 



N / N \ N / ^\ 

n dp? 6 p+ - E / n dp?^^'^ - E p? 

i=l \ i=l / i=l \ i=l ) 

/ n^^?4Q.-E^? / Vid^fHQ^-Y.^'2\ 

•'^ 1=1 \ i=l / •' i=l \ i=l I 
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X ^^,a([P?,P?])r«,([p?,P?]; [q?,q?])^QA[<li\<l?])- (3.36) 

The vertex operator f can be constructed, using, for instance, perturbation 
theory methods of quantum field theory. Phenomenological vertex operators 
can also be used. Here we consider the so-called "impulse approximation" . In 
this case 

f , = E(rr)^' (3-37) 

i=l 

(f W)^ = [G(°)]-^[G'(o)(5°))iG{o)] [G(°)]-^ . (3.38) 

Here 

(f(°)), = (2vr)^e.(p« + g»), J] ^(P^'^ ' Q^'^)[G^'\<l^'^)r\ (3-39) 

N 

and Cj is the change of z-th particle. 

Performing the integration in Eq. fl3.38p according to the definition fl3.29p . 
passing to the reference frame, where 

p+ = q+, Q±=0, /\^ = {P-Qf = -Kl = -{P^-Q^f (3.41) 

and using the transformation properties of wave functions 

<fp,a([x«,p?]) = $p^=o, Jt^^^P? - ^^""PA). (3.42) 

we obtain [5] 

(P,a I J^(0) I Q,/3) = £(P,a I J^(0) | Q,/?)^, (3.43) 

fc=i 

where, for instance, 
(P, a I J^(0) I Q, = (P+ + g+)Pfc(-A 1) 

X ^p^=,,S[^^\pf -x^^K^l^u.x^^\pf + (1 -a;('=))A^) 

x%=o,/3([^^^P?])- (3.44) 

Taking into account the normalization condition for the wave functions 
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i=i \ 1=1/"' i=i \ i=i J 

rQ. N / N \ N / N \ 

X / n d<if^ - E ^? / n dqf^^'^ Q± - E 

i=l \ i=l / 1=1 \ i=l / 

X i^M'.v'^X) '^"'gff'' I'a.Mb?.??]) = 1. (3.45) 

for A_L = we get 

N 

F(A2 = 0) = ^efc. (3.46) 

k=l 

Thus, the form factor at zero momentum transfer is normahzed to the total 
charge of the system. Note that general problems of normalization of three- 
dimensional relativistic wave functions have been considered in Ref. |83j . 

Thus Eqs. (I3.43P and f l3.44p define the form factors of a many-body system in 
terms of the light front relativistic wave functions (^{[x'^'^\'p^^]). 



4 Asymptotic Behaviour of the Pion Form Factor at Large Mo- 
mentum Transfer 



According to quark counting rules [H] , [85] asymptotic behaviour of the form- 
factor of composite system is determined by the minimal number n of elemen- 
tary constituents: 



|l-n 



(3.47) 



Well-known proof of the relation (13.470 consists in using equation for bound 
state wave function under the assumption of one meson exchange, since this 
mechanism dominates in the asymptotic region. Asymptotic behaviour of the 
wave function is then used to calculate the behaviout of form-factor at large 
momentum transfer. 



Let us write the equation for fermion-antifermion system in the one gluon 
exchange approximation: 



-•2 1 2 

p2 _ P± + rn 
xil — x) 



$p(p) 



j [dq\v{x,p±;y,q±){p^ + m)0^^'^^p{q)0^'^\-p2 + m). 



Here 



<l>p(p) = $p(x,p_l) = (pi + m)i)p{p){-p2 + 



(3.48) 
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yi^-y) 

are matrices in the corresponding fermion-gluon vertices. One gluon ex- 
change quasipotential is given above by Eq. (12.291) . 

It is known |86j that the most general form of the Bethe-Salpeter amphtude 
for pion is of the form: 



Xp{p) = 75 (Xi +PX2+PX3 + [P,P]X4), 



(3.49) 



where Xi are scalar functions of invariants Pp and . Eq. (I3.49P can be rewrit- 
ten as: 

Xp{p) = 75(6 + + M3 + PiMi)- 



Using the equation 



Pi-Pi = ^ (P- -P-)7^ 



we obtain 



Xp{p) = 75(^1 + 6Pi + -^3^2 + AP2) + 2 ^5 



(Pi,- -Pi.-)674 



2(^2,- -P2,-)UPi1+ 



+ {P2 - - P2-)^3l+^{Pl,- -Pl,-)Ul+P2 + 

Using the relation: 

(pj + m)pj = m{p^ + m) 
we get the most general form of the quasipotential wave function 

*P(P) = (Pl + m)[ipi{p) + l+V2{p)]l5{-% + "^)- 

Quasipotential equation (13.481) is now rewritten as 



P 



2 pi + 



(3.50) 



_ ^ (Pi + m)[ipi{p) + 7+¥'2(p)]75(-P2 + m) 

± / [dq\V{p,q){p^ + m)0(^^(^i + m)[^i(g) + 7+^2(9)] 
x(g2 + m)0^'^\p2 + m)75. 



"— " sign is used in the case, if 0*^^^75 = —75(9*^^^. For operators O*^*-' we use 
(9{«) = 1, 'y^, 'y^^g. In all four cases we get for the asymptotic behaviour of 
wave functions: 

~ {pir\ mp) ~ ipir'- (3.51) 

Thus, asymptotic behaviour of wave functions does not depend on the concrete 
dynamics and is of the above form. 
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Let us proceed now to the study of the asymptotic behaviour of the form- 
factor. In the impulse approximation the pion electromagnetic form-factor is 
of the form: 



{Q\J,\P) = {P + Q)^F{k^) = (2P + k)^F{e) 

[rfg]5(p2 - g2)^p{$Q(£)7i'^$p(p)} + (1 ^ 2). 



2(27r)^ 



Inserting the wave function in the form fl3.50p into this equation and integrat- 
ing over dq^ we obtain 

[dp\S,m + ^^(^l)(3^'~^ + m)7i'Ht^'^ + +7+^2). (3.52) 



Here tpi depend on pj_ + (1 — x)A_l, where 

'P 



In the reference frame, where P+ = Q^, we have: A_|_ = k± and fc^ = 
In the limit A ^ — 00 asymptotic behaviour of the form-factor is determined 
by the asymptotics of wave functions. Since asymptotic behaviour of wave 
functions does not depend on concrete dynamics, the asymptotic behaviour of 
the form-factor is model independent, as is required by the automodelity [84J, 



Calculating the "+" -component in the expression fl3.52l) . one obtains: 



.Ai)-./"^«(l-y,-,,)/d,-,. 



X 



+ 1 (g 1 + 

1/2 / 



+ (l-x)(glAx) 



iplifi + 2mP+[(x + y2)yi(p*i(p2 + xipl^pi] + Qxyiy2P+^*2^2 )■ 



If we take (^2 = (such a possibility is also considered sometimes) and use 
only (I3.5ip ip\ ~ (gj_ + (1 — x)A_|_)~^, in the limit A \ ^ 00 for the form-factor 
we get Ft, ~ (A^)~^, which contradicts the quark counting rules [Hj, [85] . 
At the same time taking into account the second structure, which behaves 
as (gj_ + (1 — a;)A_|_)~^, we get the correct asymptotic behaviour of the pion 
form-factor: 

F.(-Al)~(Ai)-\ at Al^oo. 

It is essential that this behaviour does not depend on the concrete form of 
dynamics. 
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5 Scattering of Relativistic Composite Systems 



Experimental study of high energy processes during the last decades revealed 
a number of scaling properties of observable quantities. Many of these proper- 
ties can be understood on the basis of the composite quark-parton structure 
of elementary particles. In particular, the asymptotic scaling property of dif- 
ferential cross section of hadron-hadron scattering 



da 



s— >oo S 
\t/s\=const. 



f {cos 9 s), 



(3.53) 



where N is integer number, can be explained in the framework of dimensional 
analysis and assumption on three-quark structure of baryons and quark-anti- 
quark structure of mesons (quark counting rules) [HI], [85] . 



In connection with the development of composite models of elementary par- 
ticles a problem of the description of their interactions becomes of special 
interest. Study of intersections of relativistic composite systems is important 
also in connection with current and future experiments with beams of rela- 
tivistic nuclei. Here we outline a method for the treatment of problems of that 
kind 



Below we present a description of the scattering of two composite particles. 
It will be shown that some simple assumptions on the hadron interactions in 
the scattering process allow to reproduce the results of quark counting rules. 



Consider the eight-point Green function G: 



G{xi,x2, xs, X4; yi, y2, ys, yi) 
(0 I r(v9i(xi)v?2(a;2)<^3(a;3)</'4(a;4)v'|(z/i)¥'^(l/2)<^3 (?/3)^4 (1/4)) I 0) 



[(2vr) 



Jl d!^Pi d^qi exp 



i=l 



i=l 



[{2n 



,41-8 



^4p(12) ^4^(12) ^4p(34) ^4^(34) ^4g(12) ^4^(12) ^4g(34) ^4^ 



(34) 



X exp 



^(p(12)^(12) ^ p(34)^(34) ^ p(12)^(12) ^ ^(34)^(34) 



_g(12)y(12) _ g(34)y(34) _ ^(12)^(12) 



^(34)^(34)^ 



xG'(P(12)^p(12). p(34)^p(34). g(12)^ ^(12). g(34)^ ^(34))_ 



(3.54) 



In (Km the momenta P^^^\ p^^^\ P(34)^ ^(34)^ g(i2)^ ^(12)^ g(34)^ ^(34) 
introduced according to the following equations 



d(12) I (12) Pi "~ P2 j5(34) , (34) Ps ~ Pi 

'=Pl+P2, F - ^ , '=P'i+Pi, p^ ' = — - — 



Q 



(12) 



qi + ^2, q 



(12) _ gl ~ g2 



, Q 



(34) 



93 + 94, q 



(34) 



93 - 94 



(3.55) 
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Passing now to the "two-time" description, we introduce the hght front vari- 
ables and define the quantity 

= r dp^l^^dp^^'Uq'^l'^dq'^''^G{P^''\p^'^^- P(3^),p(='^); Q(12)^ ^(12). g(34)^ ^(34))^ 

J —00 

(3.56) 

where p± = Po ±P3, Rl = (^1,^2)- 

Introduce now the quantity M by the equation 

G(p(12) ^ p(12) . p(34) ^ p(34) . g(12) ^ ^(12) . g(34) ^ ^(34) ) 

= 1 d^p^'^yd^p^'^yd^^'^yd^q^'^y 

xGu{P^''\p^''\p^'^y) Gu (P('') , P^''^ , P^'')' ) 

X M(P(12) ^ ^(12)' . p(34) ^ ^(34)' . g(12) ^ ^(12)' . g(34) ^ ^(34)' ) 

xGi2(g(^^\g(^^)',g(^^))G34(Q(^^\g(^^^g(^^)) = (Gi2G34)M(Gi2G34). (3.57) 

The quantity G can be presented in the form (here and in what follows we 
omit the arguments which are related to the relative momenta and this will 
not cause any misunderstanding): 

^(p(12)^p(34)^g(12)^g(34)) 

=Gl2(p(^^))G34(p('')).M(p(^'^g(^'^p(''^Q('')).Gl2(g(^'))G34(g('')). 

(3.58) 

The symbol • in 03.581) has to be understood in the following sense 

xi(...,pf\pf^pf\pf')B(pf',pf';pf',pf',...) (3.59) 

and dots correspond to the set of other arguments on which the operators A 
and B can depend. 

Knowing the pole singularities of the two-particle Green functions 6*12, G23 
one can show that in the vicinity of these poles the function G looks as follows: 

^(p(12)^p(34)^g(12)^g(34)) 

^ ^2.2,^4l4 ^12(P(^^))^34(P(^^))^|2(g^^^))^^4(g^^^7 

i y ) i (p(12)2 _ Mi22)(P(34)2 - M|4)(Q(12)2 _ M{1)(Q(34)2 _ M^l) 

X^+(p{12))^+(p(34)).ji^(p(12)^p(34)^g(12)^g(34)) 

.^12(g^''))^34(g(=''^). (3.60) 
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Here 



Mi234= lim Gr2^(p(^2))G3i(p(^^)) 

Q(12)2^A'/(2,Q(34)2^J\^^2 

.Gi2G3lMGi2G34(P(^^\ g(^^\ Q^''^) . Gr2^(g(^^))G3-;(g(^^)), (3.61) 

Mf2, MI4, M(|, M3I are the masses of corresponding states. 

From Eqs. fl3.58p and fl3.60p we get the following expression for the scattering 
amplitude: 

^^p(12)^p(34)^g(12)^g(34)) ^ vl>+ (p(12) ) vl>+ (p(34) ) 
.M(p(12)^ p(34)^ g(12)^ g(34)) . (Q^''^ ) ^34 (Q^'"^ ) • (3.62) 



Eq. (13.621) gives a general expression for the scattering amplitude in the case of 
scattering of composite particles. The detailed form of the scattering amplitude 
depends on the interaction mechanism in the intermediate state and on a 
special form of the wave functions of the scattered objects. 



6 Particle Exchange in the Intermediate State 



Below we consider two possible mechanisms: 1) scattering via the exchange 
of some intermediate particle (Fig. Hj), 2) constituent interchange mechanism 
(Fig. [5]). For the sake of simplicity we consider the case of scalar particles. It 
can be shown that in the first case in the reference frame, where 

p(i2) = pf ), gf ) = gf ), t = (P(^^) - g(^2))2 = -A f 



one has the following result for the scattering amplitude 



T = 



_J1 



+ Fi^ (t)Fif it) + (t)F3« [t) + F[f {t)F^f {t)] , (3.63) 
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Fig. 4. 

where -^12^ -^34^ scalar form factors of the scattered particles, which are 
expressed in terms of the light front wave functions in the following way: 



(27r)3 /-i dx 



2 Jo — x) 

X J dp^^+^,,^^^{x,p± + (1- x)Af)<!?pa2)^^{x,p±), 
{27if /-i dx 



2 Jo — x) 



X 



y(ipl$+(34)^p(x,px - (1 -a;)A?)$p(34)^Q(x,px), 



(27r)3 /-i dx 



2 Jo a;(l-a;)2 

X J o?pl$t{i2)^^(x,px -xAj^)$^(i2)^o(a;,pl), 
(27r)3 /•! dx 



2 Jo a;(l — xy 



(3.64) 



7 Constituent Interchange Mechanism 



2(27r)3 Jo 

+ (1 - x)A ?)^t^'^lix,p^)[M-t, + M 



Considering the constituent interchange mechanism (FigjS]), one gets the fol- 
lowing expression for the scattering amplitude 

T = ^[' -J^ I rfp.$t(-) (X, p. - .A i") 

34 

- S{x, p± + xAf -{1- x) A - S{x, p±)] 

X <|.(i^to(^'P"'^-^^^^)*?'Lo(^'^"'^ + (l-^)^?)' (3.65) 
where A ? = -t, A = 
Here the notation has been introduced 

5(a:,p^) = !^i±^ + ^^i±^ (3.66) 
1 — X X 
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Fiff. 5. 



and the following properties of wave functions have been used: 



= $(1 



(3.67) 



Let us consider now the wave functions of the composite particles in the form 



, N = A,B,C,D, 



(3.68) 



A, B and C, D denote the hadrons before and after the scattering and corre- 
sponding powers, respectively. 



Inserting the wave functions (13.681) into Eq. (13.651) for the scattering amphtude 
one gets in the asymptotic region 



T 



1 



\t\^oo gA+C+D-l \ 2 



1 + z\~^ fl - z 



f{z), (3.69) 



where 



dx(f\{x)ip^{x)Lpcix)LpD{x) 



:i-a;)2i^ + x2l±2 



(l_a;)2i±i + 3;2L^ 



X [x{l - a;)^+^+^+^-3x"2^(l - x) 



-2D 



(3.70) 



(2vr) = 



dp±^%{x,p±)[x{l -x)] 



-B 



z = COSIES, where -^s is the scattering angle in the c.m.s. 



s s 
-t^-{l-z), -u^-{l + z). 



Eq. (I3.69P is in close connection with the results of quark counting rules [M 
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IV. Deep Inelastic Form Factors of Composite Systems 
and Multiquark States in Nuclei 



The great interest to deep inelastic interaction processes is caused by the pos- 
sibihty of studying the internal structure of hadrons and nuclei experimentally 
and checking different theoretical models on the assumptions about compos- 
ite nature of strongly interacting particles. The main part of experimentally 
observed properties of these processes (in particular, the scale properties of 
structure functions) have been explained in the framework of composite quark- 
parton models of hadrons, in which the hadron is considered as bound state of 
some parallelly moving point-like constituents [88]-[91]. Interaction between 
constituents and their transverse motion inside hadron is neglected. 

More precise measurements in wider range of kinematic variables have led to 
the discovery of deviations from exact scale invariance in the behaviour of 
structure function [92]-|96]. There were attempts to explain these deviations 
on the kinematical (search for new scale- invariant variables [97j) and dynamic 
(taking into account higher chromodynamical corrections [98]- |100] ) basis. 

Here we incorporate the transverse motion of constituents in the composite 
system, which leads to the violation of Bjorken scaling of structure functions. 
In this approach, hadrons are considered as bound states of quarks, described 
in terms of light front wave functions. 



1 Construction of Tensor Wfj,i, 

Consider the quantity R^u, which is defined by the vacuum expectation value 
of the chronological T+-ordered product of the Heisenberg field operators 
ipi{x^^^) and local currents and J^: 



= (0 I T(^,(x«) ■ ■ -^^(xif ))J,(z)J.(0)^+(y«) ■ ■ -^^(z/i^))) 1 0) 



= (27r)-«^ / n d^Q'-'^ exp - i ^(p^^^x^^^ - q^'^y'^'^) 




xR,.{[p%[q%z). 



(4.1) 



Here [x^*^], [p'-*-'], [q^^'^] are the sets of corresponding 4- vectors. 



The quantity R^i, can be presented as (see Fig. E]) 



R, 
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ii 



R, 




Fig. 6. 

where G is the A^-particle Green function of fields Lpi{x'^^): 

G(\xf\, \yf]) = (0 I T(^i(a:«) ■ ■ ■ ^Tv(xf ■ ■ ■ ^W^^)) I 0) (4.2) 

and "two-photon" vertex r^,^ is defined by the sum of irreducible diagrams 
with 2N + 2 points (legs). 

Let us introduce now the three-dimensional quantity R^u, equating all 
Xj^ and 1/+^ = y+ in fl4.ll) 

= (0 |T(y,i(x+,xL'\fJ))---^^(x+,xL^\ffVMW'/.(0) 

M{y^J-\y'^^)---vUy^.y^^\y^^'')) I o) 



(2 



1=1 



X exp 



- 2 (^+^" ^ y+Q- 



N 

2 



i=l 



N .1 



i=l 



(4.3) 



Fourier transforms of R^^^, and Rf^^, are related to each other in the following 
way: 



R,u{[p?,P?];[q?,q^'^'- 



N / ^ \ / ^ \ 



i=l \ 1=1 / \ i=l 



Let us single out now the contribution of A^-particle bound states in matrix 
element (14. 3p . expressing the T-product via 6'-functions and using the com- 
pleteness of physical states. Using the integral representation of ^-function, we 
obtain that the quantity -R^i, has the double pole singularities in the vicinity 
of the points corresponding to bound states with masses and Mg and a 
set of other quantum numbers a and /?, respectively: 

R,.{\p?,P?];[d\<i?];z) 
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1 2 



(27r)^ 



■ (4.5) 



Here ^p([p+'',p|l'']) is the three-dimensional wave function of many-body sys- 
tem and is defined in the following way: 



I W_dx''l^ dx exp 
1=1 



N 



2 2^ I —p\.x_ — pyx Y 



X (0 I (^i(0, f • • • (^Ar(0, f I ^, «). (4.6) 



Define the three-dimensional vertex function 

^..([p?,p?];k?,gi^];^) 

„p, iV / iV \ „ iV / N \ 

1=1 \ i=i / •' i=i \ i=i I 
-rh i'ii^'s {q. - T. 9?) f n "if'/''" (qi - 1; 9?') 

^ i=l \ 1=1 / 1=1 \ i=l / 

xG(P;bJ\rf^];[pf,pf])f,.([pf ,pi)'];[gf,gf];^^ 

xG(g;[gf,gf];[g«,g«]), (4.7) 

where G is the three-dimensional Green function of iV-body system, which has 
the following pole singularity in the vicinity of bound state \P,Oi) with mass 
Mq and set of other quantum numbers a: 



1 *p([p?,P?])*$([gV^?T]) 



■ (4.8) 



Taking into account the pole singularities of three-dimensional Green functions 
at and Q'^ ^ M| we can present the quantity R^^ as follows: 



(27r)^ 



^p(K^pT])^a([gv^gT]) 

(^^-M2)(Q2_M2) 



X 



X 



•'^ 1=1 \ i=i ) ■' i=i \ %=i / 

t n (q. - f: 9?') / n •'si'''''^' [q- - z 

•'^ i=l \ 1=1 J 1=1 \ i=l ) 
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x^J([pf,J?f])f,.([pf,pi)'];[gf ,gi)'];^)^Q([gf,gf]). (4.9) 

Comparing this expression with fl4.5p . we get the following expression for the 
matrix element of T-product of currents: 

{P,a\T{J^{z)MO))\Q,P) 
-"^ 1=1 \ i=\ I i=l \ i=l / 

rQ. N / N \ N / N \ 

•'^ i=l \ i=\ ) •' i=l \ i=l ) 

x*J([pJ\pi^])f,.([pJ\p?];[gJ\g?];^)^Q([g?,g?]). (4.10) 

Fourier transform of this matrix element defines the amplitude of virtual 
Comptin scattering of photon with space-like momentum on the hadron 
with momentum P^: 

T^,(P,g) = % J d'ze''^^P,a \ T( J^(^) J,(0)) | Q,P) 
i=i \ i=i / i=i \ i=i ) 

rQ. N / N \ N f N \ 

•"^ i=i \ i=i / i=i \ i=i ) 

x^^([p?,P?]) / rfSe^^^f,.([p«,p«]; [g«,g?];^)vI/Q([g«,g»]). (4.11) 

According to optical theorem the tensor Wp,y which defines the hadronic part of 
deep inelastic lepton-hadron scattering cross section is related to the imaginary 
part of the amplitude of the zero angle virtual Compton scattering in the 
following way: 

W^,{P,q)=Y^ f dhe'''^P,a\ J,{z)J,{Q) \ Q,(3) = ^ lmT^,{P,q), (4.12) 

^ J 271 

Taking into account the current conservation, the tensor W^,y can be expressed 
via two invariant structure functions Wi and W2'- 

W,.{P,q)=(^-g,. + ^^W^iq^u) 

+ 9.) [p. - ^ g.) W,{q\ u), (4.13) 

where Mu = Pq, M is the hadron mass. 
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a h 

Fig. 7. 

Thus, using Eqs. f l4.11l) - fl4.13l) one can express the structure functions of deep 
inelastic lepton-hadron scattering in terms of relativistic many-body wave 
functions, describing the internal motion of partons inside hadron, and the 
"two-photon" vertex function: 

= / d'z / n dp?' 5 - E P?' / n dpf'S^'^ (P^-T. Pf 

^ 1=1 \ i=l / i=l \ i=l / 

xG-\P; \j^,p?]- [p?',P?'mG^G]{[p?',p?']; [gf ,gf ];^) 

xG-i(Q;[gf ,gf];[g«,g?]). (4.14) 

Here wave over the GT^^yG denotes the integration over the "— " components 
of 4-momenta p*^*-* and g*-*-* . 



2 Lowest Order in the Electromagnetic Interaction 

The "two-photon" vertex operator T^i, can be constructed using methods of 

perturbation theory and expanding the functions and R^u = GT ^^G in 
the series in coupling constant. In lowest order we obtain: 

= [G(°)]-i[G(o)r?jG(o)] {&%^ , (4.15) 

where multiplication is understood as integration over all "-I-" and " ±" compo- 
nents of 4-momenta p^*^ and g*-*-* , is the three-dimensional Green function 
of N free particles. 

In the lowest order two types of diagrams, shown in Fig. [7] will contribute in 
r|J)]. The contribution of z-th diagram a) is equal to 



flt([p?,P?];[g?,gi^];g) 

/ N \ / ^ 
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(2p« + g)^(2p« + g), 



N 



+ g- + 



— u) — 



+ 



(4.16) 



+ ie 



The contribution of diagrams of type b) is equal to: 



(«) (i) 



W_ Jfc)u(2)r^W_.-WN (2p('^ + g)M(2p(*) + g), 



X n M - dv^HpT - ) - 

A;=l g 



^ IS 

1+ 



(i) 



i=i « 



0) 



+ ie 



3=1 « 



^(»)2 , 2 ^(i)2 , 2 

(I) ^ (J) \- ^£ 

P+ Q+ 



i=l p+ 



P-- E u) ' + q~ 

j=i pX' 



pV 5+ 



(4.17) 



Assuming that the partons constituting hadron are on the mass shell and ne- 
glecting small terms of order P- — Y1, {p ± +^?) /p+\ "we obtain that diagrams 

i=l 

of type b) do not give contribution to r|°J and contribution of i-th diagram of 
type a) is equal to 



rZ{\p?,P?]^\q?,q?]^q) 



ef (2p«+g)^(2p« + g), 



^1/ _ 



2^(47r)^-ipW[(pW+g)2_^2 + 

N 

Y[p?^ip?-q?)^^'Hp?-q?)^ (4-i8) 



i=l 



where os the electric charge of z-th parton, and p is the momentum of 
parton on the mass shell: 



Summing the contributions of all diagrams of type a), inserting into the ex- 
pression (14. lip for virtual Compton scattering amplitude and extracting the 
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imaginary part of the expression obtained, we get for tensor W/j^iy-. 
1 rP+ ( ^ \ 



X 



X 



N / N \ / N \ 

Y: (2p« + g),(2p« + + - m^] . (4.19) 

i=i lp+ J 



Introducing the wave function (!?p{[x^'\p?]) (see Eq. ((233])), we can rewrite 



the tensor Wfj_i, in the form: 



4(47r)^-i Jo 



X 



„ Af / N \ 

i=i \ i=i / 



X E 1^ (2p^^^ + Q),i^p^'^ + Q)Aip^'^ + qf - ^']|- (4.20) 



Expressions for the structure functions Wi and 14^2 can be obtained multiph- 
cating the tensor W^y by the "projection" operators L^: 



(4.21) 



Using the following expressions for these "projection" operators 



p p 



M2(l - i/2/g2) 



= - 

'^^ 2(2 - I/2/g2) 



M2(l-iy2/g2) 

for the structure functions Wi and z/W^2 we obtain: 
' 8(4vr)^-iyo H V / 



X 



I. N / ^ \ 

/ n dpf5^Ap^ - Y.P? |$p([x«,p?]) 



N ( ^2 

xE 



\ 1=1 

(Mi/ + 2p«P)2 

M2(l - I/2/g2) 



(4m/ - g^) 



(4.22a) 
(4.22b) 



5(g2 + 2p«g) , (4.23a) 
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V 



K47r)^-i (1 - z/2/g2) Jo ii 3, 



1 ^ ( ^ 



.(0 



(0 



i=l 



\ i=\ 

3(Mz/ + 2pWp)2 

M2(l-Z/2/g2) 



- (4m2 - g2^ 



5(g2 + 2p»g)L (4.23b) 



For further consideration we proceed to the frame, where the virtual photon 
and hadron are moving along the z axis: 



In this frame 



P=(P_,P+,0^), g=(g_,g+,0^). 



2p«P = x»M2 + ^^^^±^ 



and (5-function in fl4.23ap . (]4.23bp can be rewritten in the form: 



5(g2 + 2p»P) = ^5 



Here we have introduced the variables = — g2 and 



9+ 



p+ M(z/ + Vi^' + ) ■ 



(4.24) 



(4.25) 



Scale properties of structure functions with respect of variables ^ are discussed 
in several papers [97], [99]- |102] . The Nachtman variables ^ is the generaliza- 
tion of usual Bjorken variable xb taking into account the hadron mass and is 
related to by the following relation: 



2x 



B 



1 + 



(4.26) 



Expressing the kinematic variable u via the variables and ^ we obtain the 
following expressions for structure functions: 



N 



\ 1=1 / 



i=i I 



2(47r)^-ie ^0 IJi x« 
/ l[dp?6^'^P^ |$p([x«,p«]) 



— m ■ 



Pi 



(4.27a) 
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X 



i=l 

i=l I 



X(5 



^(j)2 , 2 



(4.27b) 



If one neglects the masses and transverse momenta of partons (m? <^ Q^, 
p±^'^ <^ Q^), the (5-function takes the form 



Then the structure function Wi vanishes and for the structure function W2 
we obtain 

X / n^P?^^'M Ep? E^M^^'^-O- (4.28) 



In the asymptotic limit [y^Q^ ^ M^, 0;^ is fixed) the variable ^ coincides 
with the Bjorken variable xb and we obtain that in this limit the structure 
function VW2 is scale-invariant with respect to the variable xb'- 



N 



— x^) 



i=l 



N / N \ „ N 

X / UdP±^^"HEP?]\M[x^'^:Pf])\ E^lSix^'^-XB). (4.29) 

i=l \j=l / i=l 



Normalization condition for wave function ^([x^'^\pj^]) is of the form 



1 ^ 



N 







•^u i=l -f^^ \ i=l I i=\ \ i=\ 



X 



i=l 
N 



.. N 



N 



1=1 
N 



1=1 

N 



,., gG-(P;lpg.Pa;bf.pf]) ,) = 1. (4.30) 



i=l 



9P2 
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Assuming that the interaction kernel does not depend on the total energy and 
using explicit expression for the Green function of free particles we obtain 
the following sum rule: 

— dxB = l^ei. 4.31 

Mxb ~r 



3 Model Parametrizations of Wave Functions 



Let us consider now the case, when the hadron consists of two constituents. 
This case corresponds to meson, which consists of quark and antiquark. We 
will neglect contributions of gluons and quark- antiquark sea. 

Expressions for structure functions in the case N = 2 have the form: 



STrf Ja Xil — X) J 



X 



m 



8< 
Q^x{x~i) 

e 



'^2 , 2 Q'^X{X-S) 

P± + m — 



X 



X 



8< ^ (Qye + M^y 



6Q'^x{x - 



+ Q2 _ 



P± + rn — 



(4.32a) 



(4.32b) 



Here we assume that masses of constituents are equal to each other: mi 

7712 = fTT'- 



In fl4.32ap . fl4.32b|) the limit of integration over x is defined from the 5-function 
and equals to: 



(4.33) 



Neglecting masses and transverse momenta of quarks we obtain that the struc- 
ture function Wi vanishes and the structure function takes the following 
form: 



(4.34) 
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If we choose the following parametrization for wave function $p 

^p{x,p±) = C 



'p \ + m? 
xil — x) 



a 



then for the structure function VW2 we get 



X 



8 {Qye + M2)2 

(2n-l)[m2-a^(l-^)]2"-i ' 



(4.35) 



(4.36) 



In the Bjorken limit (Q^ ^ 77^2^ ^ x^) for n = 1 we obtain |103j : 



— axB(l — xb) 



(4.37) 



Let us consider now another parametrization of the wave function 



$p(x,pj_) = C exp 



x(l — x) 



(4.38) 



Inserting this wave function into fl4.32al) and (]4.32bp and neglecting quark 
masses [jr? -C Q^), for the structure functions we get: 



2 ve? + e^)g2(l-0|^|2 



yW2{Q\i) 



l + 2/3:f^jexp(2/3^)E^( -2/3^^ 
6Mg2 g2/^2 _ ^2 



+ 



e (gve + M2)2 

(e? + e2)MQ4 Q2/^2 _ ^2 



4^2 



(Q2/^2 + M2)2 



ICI 



x|-exp^2,|,Ji..^-2^|,^ 
Here is the integral exponent function. 



(4.39a) 



(4.39b) 



Note that more realistic consideration needs inclusion of spin and color de- 
grees of freedom of constituents, which can be done directly by means of 
corresponding formalism, developed in the previous sections. 
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4 Quark Degrees of Freedom in Nuclei and the EMC-EfFect 

As has been mentioned in the Introduction, recent investigations have led 
to the conclusion that the consideration of nuclei as the systems of quasi- 
independent nonrelativistic nucleons is incomplete and they require a rela- 
tivistic description of the internal motion of nucleons in nuclei and taking 
into account quark degrees of freedom. These are first of all the prediction 
1 104] and observation |105] - [T07] of the cumulative production of particles in 
hadron-nucleus and nucleus-nucleus collisions and change of an exponential 
fall-off of form-factors of light nuclei at small momentum transfers to a power- 
law fall-off at large momentum transfers |108j . |109] according to the quark 
counting rules [H], [85] . 

In this sense the most impressive are the results of deep inelastic lepton- 
nucleus scattering experiments |110j - [TT^ (the so-called EMC-effect). Differ- 
ent models are suggested for the explanation of this phenomenon (see, e.g., 
Refs. ^1T5] - P^y] ). It seems that the above mentioned regularities are of the 
same nature and are determined by the possibility of formation of multiquark 
configurations (multiquark bags) in nuclei. Here we analyze the EMC-efFect 
and show that the effect can be explained by taking into account the scatter- 
ing on colorless multiquark configurations which are contained in the medium 
of spectator nucleons. 

Let us consider deep inelastic scattering of charged lepton on a nucleus A. 
We shall assume that in the nucleus, together with nucleons (three-quark 
bags), the configurations with six, nine, etc. quarks are formed with definite 
probabilities and leptons interact with the nucleus by the exchange of virtual 
photons with quarks from these bags. Then the nucleus structure function can 
be represented by the sum: 

F2{x, Q') = J2 NiA, K)F^{x, Q'), (4.40) 

K=l 

where is the structure function of the nucleus A which contains a 3K- 
quark bag and {A — K) nucleons. The coefficients N{A, K) of these structure 
function have the meaning of the effective number of SiiT-quark bags in the 
nucleus A and obey the normalization condition: 

A 

Y,KN{A,K)=A. (4.41) 

K=l 

We use the parametrization of N{A, K) in the form of the Bernoulli distribu- 
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tion: 

N{A,K) = (4-42) 

For the parameter p{A) determining the probabihty of a three-quark nucleon 
to get a Si^-quark bag we consider two possibihties: 

1) p{A) is determined by the ratio of the bag and nucleus volumes: 

j)(A) = YI£ = !]l = !i (4 43) 



Taking the bag radius rx close to the nucleon radius rj^- 0, 8 fm and Rq = 
1,4 fm tlSOj we obtain p{A) = 1, 1867 

2) p{A) is determined by the ratio of the bag and nucleus cross sections: 

2 

p(A) r-^^r^ (4.44) 

Ra 



The coefficient of proportionality in (14.441) for this case has been obtained in 
Ref. [131j by fitting the data on production of vr-mesons with large transverse 
momenta in photon-nucleus scattering, and we shall also use this parametriza- 
tionp(A) = 0,085A-0'57. 

In both cases N{A, K) are rapidly decreasing functions of K and the main 
contribution to the structure function is given by first few terms of the sum 

We proceed now to calculation of structure functions . As it is well-known, 
the structure function F2 appears in the decomposition of the deep inelas- 
tic tensor W^^ into gauge-invariant structures. The tensor W^i, itself is pro- 
portional to the imaginary part of the forward virtual Compton scattering 
amplitude: 

T^,{PA.q) =ijd^z exp{tqz){PA \ T{J^{z)MO)} I Pa) (4.45) 

that can be expressed through the two-photon vertex function and relativistic 
wave functions of composite systems: 

T^,{PA.q)=i r^'^ n dpfdPK,+ f n dpfdPK,± 

"'0 i=K+l i=K+l 

4pa,^-Pk,^- E p^f\5^'^{pA,^-PK,^- E p?) 

\ i=K+l ) \ i=K+\ ) 



X 
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X ['^^'^ n dqfdQK,+ I n dqfdQK,± 

•'0 i=K+l ■' i=K+l 

V i=K-^l J \ i=K+l 

x'^QAQK.^QK.^AQf.^f])- (4-46) 

Here T is the two=particle vertex functions, \E'p^ is the relativistic wave 
function of the composite system consisting of a SfT-quark bag and {A — f^) 
nucleons, Pa,^ (a* = 0,1,2,3) is the 4-momentum of nucleus A, Pk,ij, and 
Qx.iU are the 4- momenta of the SX-quark bag, p^^\ are the 4- momenta of 
nucleons. The square brackets in the arguments of and F^,^ denote the 
sets of corresponding variables: 

b?,P?] = pr^pr^^;pr^pr^^; • • • j:\p'^'- (4.47) 



Representing now the relativistic wave function ^p^{Pk,+, Pk,±, \p+ ,P±]) in 
the form 

. 3K / 3K \ 

X J iidkj,+si^PK,^-j:k,,+j 



. 3K / 3K \ 

X / n dk.^^S^'^ Pk,± - E h± *3/.([%,+, (4-48) 
i=i V j=i ' 



where and kj^i are components of the quark momentum in the SX-quark 
bag, and calculating the function F^j, in the lowest order in electromagnetic 
interaction, we get the following expression for the structure function F^: 



rl Xi ( ^ 

•'O i=K+l -^^ V i=/s:+i 

/ n dpi,±dPK,±S^^^(pK,±+ E 

i=K+l \ i=K+l / 

„ 3K / \ 

X / Udk±S'^'HPK,±-T.h±)\'^3K{[zi,ki,±])\' 
•' 1=1 \ 1=1 J 



52 



+ 2m] + 



[U/Zk] 



Zj{Zj-U/ZK) 



X 



5(2) 



2 Q'^Zj{Zj — ^A/Zk) 



{U/Zk) 



(4.49) 



Here = —Q^ is the 4-momentum transfer squared. The variable is defined 

as 



6 



2XA 



i + (i + 4Mixyg2)i/2 ^ 

where xa = Q"^ /2MaI'i v = (Paq) /Ma, Ma is the nucleus mass. The variable 
xa varies in the interval < < 1 and is related to the Bjorken variable 
X = Q"^ /2Mu by xa = {M/Ma)x (M is the nucleon mass). It is evident that 
< X < Ma/M ^ A and 



(4.50) 



Ma a ' ^ 



2x 



l + (l + 4M2a;Vg2)i/2 ' 



(4.51) 



In (14.491) Cj and rrij are the electric charge and the mass of an j-th quark, 
respectively. 

The wave functions ^A-Ki\xi,Pi,i}) and (^2,Ki\zj, A;j,±]) are related to the func- 
tions a-k{^ i\) and \&3x([fcj]) by the following formulas: 

^A-K{[^^,P^A) = {PA,+)^-''{ H X^^A-K{^i]) (4.52a) 

\i=K+l I 



The variables Xi and zj are defined as: 



Xi = , < Xi < 1, Xi = l- Zk, 

i=K+l 

h. 



The variable Zk is the ratio of "+" -components of the 4-momenta of the 3K- 
quark bag and nucleus A, pi^x and kj^_[_ are transverse momenta of nucleons 
in the nuclei and of a quarks in the bag, respectively. 

Let us choose now the wave functions $A-x([a^i,Pi,±]) and ^^K{\zj,kj^x\) in 
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the following form: 

^A-K{[xi,Pi,±]) ~ exp 

— * 

*3ii-(h-,%±]) ~ exp 



i=K+l 



3K 1.2 



(4.53) 
(4.54) 



Mi are the nucleon masses, m,- are the quark masses. 



In the deep inelastic limit {Q ^ , rrij, xa] ^ x) the Q -dependence 

in the structure functions F.^ disappears, and after the corresponding calcu- 
lations from (14.491) for the quark distribution function f^{x), which is defined 
by the following relation 



F,^{x) = {epxf^{x), 
({Cg) is the average value of quark charges squared), we obtain [127] : 



3K Ik{xa) 



A J^dxAlxixA) 



Xa 



X 

A 



K = 1,2,...,A-1, 



(4.55) 



(4.56) 



where 



IxixA) 



dZ 



1 - Zf^^-^-^ 



Ixa"' ^ I + {PK/aA){l - Zk 

a;^(l-a;A)(^-^)+('^"') 
1 + (/3x/aA)(l - Xa) 



1 



Xa 



3K-2 



Zk , 

B{A-K, 3K-1] 



X 2Fi\3K-1,A-K,1,A + 2K-1;1-xa: 



il + pK/aA){l-XAy 
1 + {(3k/(Xa){1 - xa)_ 



(4.57) 



Here B is the Euler beta-function, 2-^1 is the hypergeometric function of two 
variables (the Appel function). 

A characteristic feature of the structure function F^ is possible existence 
of a superfast quark in a nucleus which in expreme situation takes all the 
momentum of nucleus. 



The structure functions (I4.55I) - (I4.57I) are normalized to the number of quarks 
in the corresponding bag 



rA 

/ dxf^{x) = 3K. 
Jo 



(4.58) 



One can relate the parameters aA and f3K of the relativistic wave functions 
^A-K{[xi,Pi,±]) and $3i<-([zj, kj^±]) to the radii of the nucleus A and SX-quark 
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bag, respectively, 



R\ ~ 8Aa^, (4.59a) 
r\ ~ 2AK(3k. (4.59b) 



As can be seen from ( 14.561) and fl4.57p the structure functions depend on 
the ratio 

(3k a rl rlA^/^ 



aA SK R\ mlK ■ 
Choosing = 0.8 fm, i?o = 1, 2 fm we get (3k /aA = 0.109 A^^^/K. 



(4.60) 



In the experiments |110j - [m] the ratio of structure functions F2{Fe)/F2{D) is 
measured in the region x < 1 that corresponds to the region of small xa < 1/A. 
It is evident from f l4.55l) -( l4.56l) and it can be checked by straightforward nu- 
merical calculations that the structure functions weakly depend on the param- 
eter (3k/ 01 A in this region. That is why in the numerical calculations we have 
supposed that the ratio (3k /o- a is the same for all bags. For the iron nucleus 
we have taken {(3k/ aA)Fe = 0.1. But for the deuteron (considering that the 
average distance between nucleons in the deuteron is somewhat larger than in 
other nuclei) we have chosen the value {(3K/oiA)D = 0.05. 

Experimental data on the ratio F2{Fe)/ F2{D) and the curves calculated by 
the formulas fjOOj) . fj4l2|) . fH755l) - fH37l) with two parametrizations for P{A) 
are shown in Fig. [HI The curves reproduce rather well the experimental data 
in the region x > 0.2 and somewhat differ from the data in the region x < 0.2. 
This difference is probably due to the fact that we restrict our consideration 
by the valence quarks and do not take into account contributions of the sea 
quarks and gluons. 

The analysis performed shows that deviation of the ratio Fi{Fe)/F2{D) from 
unity is due to a larger contribution of the multiquark configurations in the 
iron nucleus as compared to the deuteron. Note that according to fl4.42p the 
probabilities of formation of six-, nine-, etc., quark configurations in heavy 
nuclei are larger than in light nuclei. In the deuteron the contribution of a 
six-quark state is only a small admixture to the contribution of a two-nucleon 
state. The deviation from unity of the ratio of the sum of two-nucleon and 
six-quark contributions to the contribution of a pure two-nucleon state does 
not exceed 5% for the deuteron in the whole region < a; < 2 (see in this 
connection also |132] . [133] ). 



In Refs. [21], [25] the possibility was established for extracting information 
on the quark parton functions of nuclei from the data on cumulative pion 
production, and in Ref [134] similarity in the x-behaviour of the ratio of pion 
production cross sections on different nuclei and the ratio of the deep inelastic 
structure functions of the same nuclei was pointed out. It should be noticed 
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Fig. 8. The ratio of the structure functions of iron and deuteron. Curves 1 and 2 - 
calculations with parametrizations p{A) A'^^ and p{A) ~ A~'^^^, respectively. 
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Fig. 9. The ratios of structure functions of different nuclei in the whole region of 
variable x {0 < x < A) 

that the data on cumulative production allow us to investigate the region x > 1 
(which is not yet reached in the experiments on deep inelastic lepton-nucleus 
scattering). We have calculated the ratios F2{Fe)/ F2{D), F2{Fe)/ F2{He) and 
F2{Fe)/F2{Al) in the whole region of variable x {0 < x < A). The curves of 
these calculations in the parametrizationp(y4) ~ A~^ in the double logarithmic 
presentation are given in Fig. [91 The similarity in the behaviour of these curves 
and the ratios of cross sections of cumulative pion production is observed. 

The analysis performed shows that the explanation of the observed deviation 
from unity of the ration Fi{Fe)/ F2{D) requires the consideration of mul- 
tiquark configurations in nuclei. It seems very interesting to analyze the A- 
dependence of the ratio F2{A)/D2{D) contained in flOOD . fICTD . flT5D - (l437l) 
(there are the data on this dependence [113j ) and carry out experimental in- 
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Fig. 10. The ratios of structure functions of '^He, ^^C, '^^Ca, ^^^Ag nuclei to struc- 
ture function of deuteron versus variable x. Data from Ref. |113j . 




Fig. 11. The A-dependence of structure function ratios at various values of variable 
Q2. Data from Ref. [TT3] . 

vestigations of the deep inelastic lepton-nucleus scattering in the region, where 
the variable x > 1. It seems that until now in the region 1 < x < 1, 4 only the 
structure function of ^^C nucleus has been measured [135]. 

On Fig. [TD] the theoretical curves on the ratio Fi{A)/F2{D) for several dif- 
ferent nuclei and experimental data on deep inelastic scattering cross section 
ration aA/crr) are presented and good agreement is observed. More correct 
comparison of these quantities needs taking into account the contributions 
from structure function Fi to the cross section. On Fig. [11] the A-dependence 
of these ratios is presented at various values of variable Q^. The Fig. [T^ shows 
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Fig. 12. The 74-dependence of ratio F2{A) / F2{Pb) at different values of variables x. 

the A-dependence of ration F2{A) / F2{Pb) at different values of x. Quantitevely 
these curves at x > 1 have the same behaviour as the data on cumulative pion 
production. 

Note that in the region of very small x the experimental situation has been 
changed significantly |136j . |137j and the problem of the theoretical expla- 
nation of this phenomenon is of special interest (see in this connection Ref. 



V. Processes Involving High Energy Nuclei and Problem 
of Relativisation of Nuclear Wave Functions 

1 Scale-Invariant Parametrization of the Deuteron Relativistic 
Wave Function 



In this Chapter we consider the relativization of nuclear wave functions. Let 
us start our consideration with simplest case of a two-nucleon nucleus, i.e. 
deuteron. To construct a relativistic wave function of deuteron, we use the 
fact that in the framework of light front formalism the Lorenz- invariant com- 
bination + m'^)/x{l — x) (where p± is the relative transverse momentum 
of the constituent, x = 1/2 + (po + Pz)/{Po + Pz), Pi^i and P^, are the rela- 
tive 4-momentum of the internal motion of the constituents and the total 4- 
momentum of the composite system, respectively, m is the constituent mass) 
plays the role analogous to that of p"^ - the square of three-dimensional rel- 
ative momentum. The relativistic wave functions can be obtained from the 
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corresponding nonrelativistic expressions by substitution [6], |139] . |140] 

a;(l — X) 

and changing the nonrelativistic numerical parameters by the relativistic ones. 

In this way one can obtain, for example, the relativistic analog of the Hulten 
wave function: 

(5.1) 

which is written in an arbitrary reference frame for the arbitrary momenta 
of the deuteron as a whole and arbitrary intrinsic momenta of its constituent 
nucleons. Note that relativization of other, more accurate, deuteron wave func- 
tions does not meet any principal difficulties. In particular, in what follows we 
shall use the relativistic analog of the Gartenhaus-Moravchik wave function. 
We shall not describe here other possible ways of relativization of the wave 
functions, which can be found in Refs. |141j - [T7r] . 

In the deuteron rest frame, when the momenta of internal motion of its con- 
stituent obey the condition \p\/mM <^ 1, the wave function (15.11) turns to the 
well-known nonrelativistic Hulten wave function 

^nr{p) = CNRip' + a%j,r\p' + PIj,)-\ (5.2) 

In f l5.ip and fl5.2p a^, (3r and ajqR, Pnr are adjustable parameters of the 
relativistic and nonrelativistic wave functions, respectively, Cr and Cnr are 
normalization coefficients. 

Neglecting the dependence of variable x on binding energy from the condition 
that in the nonrelativistic limit the wave function goes over into the wave 
function ^nr we obtain the following relation between the parameters 

aR = A{m^ - alji); /3r = Aim^ - plj^). (5.3) 

From the normalization condition for wave function (15.21) 

Jdp\^^Rip)\' = l (5.4) 
we get the following expression for normalization coefficient C^r 

Cnr = -{aNR + pNRY^'^aNRpNw (5.5) 

71 
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^Rix,p±) = Cr 



xil — x) 



OiR 



xil — x) 



-A 



R 



For the normalization of the wave function (15.11) it is necessary to know the 
form of all interaction within the two-particle bound state. Assuming that 
the total interaction kernel does not depend on the total 4-momentum of the 
deuteron, we obtain the following normalization condition: 



Inserting into this normalization condition the wave function (15. ip we get the 
following expression for normalization coefficient Cr 



Cr = 2^/'{aR - Pr) \f{aR, Pr) + fiPR, aR) 



-1/2 



(5.7) 



where 



f(,aR,pR) 



4:[m'^{aR - (3r) + afl(4m^ - or)] 
iaR-PR)af{Am^-aRy/^ 

X arctg 



— aR 



1 



(5.^ 



The relation between the wave functions ^r and ^nr in the nonrelativistic 
limit with account of (15.31) takes the form 

$r(x,p±) ^27ri/W/2$^^(p). 

Here, the relativistic wave function ^r is normalized by the condition (15.61) and 
the nonrelativistic wave function ^nr is normalized by the condition (15. 4p . 



2 Break-up of the Relativistic Deuteron and the Verification of 
Scaling Properties of Its Wave Function 



Let us consider the interaction process of relativistic deuteron with target in 
which incoming deuteron breaks-up and a system of hadrons is created. 
In the impulse approximation this process is described by two diagrams in 
which only the one nucleon of the deuteron interacts with a target. Another 
nucleon, the so-called spectator-nucleon, does not interact with the target and 
observation of its distributions allows obtain information on the character of 
internal motion of constituent in deuteron to be obtained. In the range of 
small transverse momenta this approximation seems very reasonable. 

For the study of the dynamics of the relativistic deuteron experiments with 
hydrogen targets are the most convenient, since in this case there are no effects 
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associated with disintegration of the target, and the selection of the spectator 
nucleons makes it possible to obtain direct information about the deuteron 
wave function. 



(5.9) 



We consider now the process of interaction of the relativistic deuteron with 
hydrogen target in which a spectator-nucleon and the system of hadrons 
are produced Z) + p — >■ Pspi^sp) + Xjsf. One of the simplest processes of this 
type is the process of direct break-up of the deuteron Dp — > ppn. 

Assuming that the spectator-nucleon does not interact with the target we 
can obtain the following invariant distribution for spectator-nucleon |6], [140j . 

da Ai/2(sAr7v,m2,m2) ^r{x,p±] 

Here s is the usual Mandelstam variable for the deuteron-proton system, snn 
is the similar variable for the subsystem consisting of the interacting nucleon 
of the deuteron and the target proton. Energy-momentum conservation leads 
to the following relation between these variables: 

Snn = s{l - Xsp) +m^- ^^P'^^"^ , (5.10) 

(^in{sNN) is the total inelastic cross section of the nucleon-nucleon interaction 
in the given channel A^A^ X^, rn is the nucleon mass, m/) is the deuteron 
mass. This flux factor is defined in the following way: \{x,y,z) = {x — y — 
zY — 4yz. The variable Xgp is defined as follows: 

Y — {Esp + Psp,z) (^'\'\\ 

{E, + P,,.) + {Ep + Pp,,)- ^"-''^ 

Here Egp, Ed, Ep and Pgp^z, Pd,z, Pp,z are the energies and the 2;-components 
of momenta of the spectator nucleon, colliding deuteron and proton, respec- 
tively. Note that the variable Xgp is scale-invariant and Lorenz-invariant under 
transformations of the reference frames along the collision axes (the z-axes). 

The arguments of the wave function $_r(x, p±) are related to the variables X^p 
and Psp,± in the following way: 

^ = l-fl + f^^)^.p, (5.12) 

\ + ^d,z / 

pi = — Psp,± in the reference frame, in which Pd^±^ = . 



It follows from these expressions that the experimental observation of the 
spectator-nucleon distributions allows information on the internal motion of 
nucleons inside the relativistic deuteron to be obtained. 



61 



In the high energy hmit the distribution of the spectator, summed over all 
possible hadron systems X^v, takes the form 



1 da 



\^r{Xsp, Psp,±) 
1 — Xsp 



(5.13) 



which is very close to the predictions for inclusive distributions obtained in the 
framework of the hypothesis of limiting fragmentation ^149] and those of the 
parton model [88]-[9T] and of the automodelity principle for strong interaction 
|150] . [151j . Deviations from the automodel behaviour may occur because of a 
possible weak dependence of wave function parameters on the energy. So, the 
experimental study of processes with beams of deuterons of different energies 
seems very interesting. We shall discuss this point later on when comparing 
theoretical calculations with experimental data. 

A characteristic feature of the spectator-nucleon distribution (15.91) in the ref- 
erence frame, in which the target proton is at rest, is the prediction of a 
maximum at the value of variable Xgp 



Ed + Pd,: 

which approaches its limiting value Xgp = 0,5 as energy increases. As it will 
be seen from what follows the positions of the maxima of the experimental 
X^p-distributions agree with the values predicted by (15.140 . 

To compare the relativistic parametrization (15.11) with the nonrelativistic Hul- 
ten wave function (15.21) . we consider the momentum distribution of the spec- 
tator nucleons in the deuteron rest frame [152] . The momentum distribution 
of the spectators is related to the invariant differential cross section (15. 9p in 
the following way 



da 

dP7p ^ {PI + m 




(5.15) 



In Fig.[T3]the results of theoretical calculations with the relativistic wave func- 
tion (15. ip and with the nonrelativistic Hulten wave function (15. 2p are compared 
with the experimental distribution of spectator neutrons in the process of di- 
rect break-up of the deuteron D+p^p + p + Ugp in the deuteron bubble 
chamber bombarded by deuterons with 3,3 GeV/c momentum and converted 
to the antilaboratory frame in which the momentum of the incident protons 
is equal to 1,66 GeV/c jT53]-[T55]. 

Experimentally, the spectators were chosen as particles having the smallest 
momentum in the deuteron rest frame. The small admixture of other nucleons 
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Fig. 13. Momentum distribution of the spectator neutrons in the deuteron rest 
frame in the reaction Dp ppngp. The solid and dashed curves are the results of 
calculation with the relativistic and the nonrelativistic wave functions, respectively. 

cannot affect significantly our considerations. Numerical values of the param- 
eters aR and Pr of the relativistic wave function of deuteron were calculated 
by Eq. (E3]) from the values onr = 0,0456 GeV/c and Pnr = 0,26 GeV/c 
|156] of the nonrelativistic Hulten wave function. The cross section of elastic 
nucleon-nucleon scattering at these energies is almost constant and equal to 
24 mb [157] . [158] . The solid curve on Fig. [T3l corresponds to the theoretical cal- 
culation with the relativistic wave function (15. ip and the values of the param- 
eters an = 3,521 {GeV/cf and (3r = 3,390 (GeV/c)^ {x^/Np = 27/20). The 
dashed curve corresponds to the calculation with the nonrelativistic wave func- 
tion (15.21) and the parameters unr = 0,0456 GeV/c and Pnr = 0.26 GeV/c 
{x^/Np = 65/20). Comparison with experimental data shows that in the range 
of spectator momenta Pgp < 0,2 GeV/c the relativistic wave function gives 
somewhat better description of the experimental data than the nonrelativistic 
one. 

In Fig. [Hlthe momentum distribution of the spectator protons in the reaction 
D + p p + (pn^) + Psp at the energy y/s = 52 GeV is presented in the 
deuteron rest frame. The data were obtained in an experiment at the CERN 
ISR with the colliding deuteron-proton beams [36], [37]. In Fig. [15] the mo- 
mentum distribution of the spectator protons in the same reaction and at the 
same energy is presented in the frame of the colliding beams. The curves in 
Figs, [m and [15] correspond to the calculations with the wave function (15. ip . 
The cross section ain{sNN) of neutron diffractive dissosiation n+p (p7r~)+p 
at the considered energies was taken to be constant and equal to 185 mb [36]. 
The numerical values of the parameters and (3r obtained by fitting the ex- 
perimental data and the corresponding values of /Np are given in Table 1. 
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Fig. 14. Momentum distribution of the spectator protons in the deuteron rest frame 
in the reaction Dp — > p{pT^~)psp- 



500 




Fig. 15. Momentum distribution of the spectator protons in the frame of the colhding 
beams in the reaction Dp — > p{p'K~)psp- 

The agreement of the results of theoretical calculations with the experimental 
data confirms the validity of using the relativistic parametrization ( 15. ip . 

In order to verify the scaling properties of the relativistic wave functions, 
we have compared the theoretical calculations with the experimental distri- 
butions of the spectator nucleons at various momenta of incident deuterons 
[159] . The distributions of the spectator neutron calculated using the wave 
function (15. ip were compared with the experimental distributions of neutron 
from deuteron stripping, which were obtained in the following way. A beam 
of deuteron accelerated at the JINR Synchrophasotron to momentum Po col- 
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Table 1. Parameters of the Relativistic Hulthen Wave Function 



Reaction 


a^, (GeV/c)2 


(GeV/c)2 


xVNp 


Calculation by Eq. (j5.3p 


3,521 


3,390 




D+p^p+{p7r-)+Psp 








^ = 52 GeV 


3, 522 ±0,006 


3, 383 ± 0, 033 


11/32 


Deuteron rest frame 








D +p ^ p + (p-K-) + Psp 








= 52 GeV 


3, 513 ±0,005 


3,478 ±0,019 


28/14 


Center-of-mass frame 








D + Al^nsp + Xn 








Pd = 3, 46 GeV/c 


3,5156-fixed 


3, 5240 ± 0, 0004 


1,9 


Pd = 4, 46 GeV/c 


3,5156-fixed 


3, 5222 ±0,0006 


1,1 


Pd = 7, 66 GeV/c 


3,5156-fixed 


3, 5203 ±0,0010 


2,9 


Pd = 10 GeV/c 


3,5156-fixed 


3, 5170 ±0,0021 


1,1 


D + p ^ p + p + Usp 








Pd = 3,3 GeV/c 


3,5156-fixed 


3, 4572 ± 0, 0005 


232/150 


do" / dXsp 








in three intervals of Psp,± 








D + p ^ p + p + Usp 








Pd = 3,3 GeV/c 


3,5156-fixed 


3, 4579 ± 0, 0005 


311/150 



da/dPsp^± 
in three intervals of X, 



lided with an internal Alluminum target. Neutrons scattered from the target 
under the angles 6 < 0, 001 rad were directed in the 100-cm bubble chamber. 
This set-up gives with a good accuracy the neutron with zero transverse mo- 
mentum P<jp,±- Momentum spectra of neutrons were measured in the reaction 
np — > ppn~, which can be identified clearly in the bubble chamber. Events 
selected by this method were used for obtaining the momentum distributions 
of neutrons detected in the 100-cm bubble chamber. Detailed description of 
the experimental procedure can be found in [160j - [T62] . 

Assuming that deuteron break-up process occurs on the separate nucleons of 
target nuclei one can suppose that distribution of spectator nucleons does not 
depend on the target. Independence of the spectator fragment distributions on 
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the target for this type of processes seems to be general property and allows 
one to use theoretical distribution fl5.9p in this case. 



In Figs. [16] a,b,c,d the experimental distributions 



da 



sp(l\ Lgp p 



of spectator neutrons for four different values of the incident deuteron energy 
and the corresponding theoretical distributions normalized to the unity in the 
points of their maxima are presented. The relation between this distribution 
and the invariant distribution (15. 9p at zero value of transverse momentum of 
the spectator has the form: 



Performed analysis shows that in the considered region of momenta of the 
incident deuteron the parameters of the relativistic wave function and (3r 
depend weakly on the energy of the incident beam (see Table [2]). This fact 
indicates that in wide range of energy the relativistic wave function of deuteron 
$ij contains no other dependence on the energy apart from the dependence 
on the variable x and allows to say that wave function with the scale-invariant 
parametrization of the "longitudinal motion" of constituents in terms of light 
front variables gives a good description of the relativistic deuteron. 

The positions of maxima of the spectator neutron distributions in Figs. 
a,b,c,d coincide rather well with the values of Xsp predicted by Eq. (I5.14p . 

In order to study the transverse momentum distributions of the spectators 
|148j we have used the experimental data on the interaction of deuteron with 
momentum 3,3 GeV/c in hydrogen bubble chamber. From the experimental 
point of view, quite comprehensive data are available for the direct break-up 
channel dp ^ ppn. In this case ain{s]yN) in Eq. (15.91) must be replaced by 
the total elastic cross section aei{sNN) of the nucleon-nucleon interaction. We 
restricted ourselves to the case of the neutron-spectator, since there are much 
more experimental data on a^f than on a";^ (see, e.g., |157] ). It should be noted 
that in the experiment with a beam of 3,3 GeV/c deuterons, Spp varies in an 
interval, in which aei{spp) is almost constant and equal to 24 mb. 

The experimental distributions da/dX^p and da/dPsp^± in the rest frame of 
the target proton were analyzed. They are related to the invariant differential 
cross section as follows: 





(5.16a) 
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Fig. 16. The distributions 

eut( 
da 



da 



dyC spd^ sip 

menta of the incident deuterons. 

rX 



sp,. 



of spectator neutrons at different mo 
da 



da 



dX gpdPgp,. 



27r 



sp,. 



da 



Xsp J dPsp/E. 



(5.16b) 
(5.16c) 



sp 



The experimental distributions da/dXsp {da/dPsp,±) integrated in three dif- 
ferent intervals of Psp,± i^sp) were compared with the results of the theoretical 
scheme described above. The results of the analysis are given in Figs. [TTla.b.c 
and Figs. [TS] a,b,c. The theoretical curves in these Figures correspond to the 
values of parameters of the relativistic wave function (15.11) given in Table 1. 
Values of /Np (where Np is the number of experimental points) allow to 
consider the agreement of the model with the experimental data to be sat- 
isfactory. In accordance with the prediction (15.141) . the maximum in the Xgp 
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Fig. 17. The da/dXsp distributions of spectator neutron in the reaction Dp ppugp 
in the intervals a) 0.01 < Psp,± < 0.04 GeV/c, b) 0.04 < P,p,_L < 0, 07 GeV/c, 
c) 0.07 < Psp,± < 0.1 GeV/c. 

distribution is at X,„ ^ 0.44. 



3 The Break-up of More Complicated Nuclei 

It is interesting to attempt to find relativistic analogues of the wave functions 
of more complicated nuclei and to verify the universality of scale-invariant 
properties of these relativistic wave functions in the framework of light front 
formalism [S], [H], |163j - [TB5] by comparison of corresponding results with 
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Fig. 18. The da/dPsp,± distributions of spectator neutron in the reaction 
Dp ppusp in the intervals a) 0.40 < Xgp < 0.43, b) 0.43 < Xsp < 0,46, 
c) 0.45 < Xsp,± < 0.49. 

experimental data. 

In this section the process of knocking out the nucleon from the relativistic 
nucleus in the collision with hydrogen target (proton) is considered. The rel- 
ativistic nuclear wave functions having scale-invariant properties and turning 
in the nonrelativistic limit to the well-known Gaussian nuclear wave func- 
tions are introduced. The invariant differential cross section of the nucleon 
knock-out reaction from nucleus is expressed via the overlap integral of these 
wave functions. This overlap integral is calculated in an analytical form using 
the Gaussian parametrization of the light front wave functions. The results 
of theoretical calculation are compared with experimental data on the partial 
break-up of '^He nuclei with momentum 8,56 GeV/c and ^He nuclei with mo- 
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mentum 13,5 GeV/c in collision with proton target. The data obtained on the 
100-cm bubble chamber of the JINR [166], [167] and the data on the process 
^He+p ^He+X at different value of momentum of incident nuclei obtained 
at SATURNE [M] have been used. 

Let us consider the process of the knocking out one nucleon from the relativis- 
tic nucleus A on the hydrogen target. Assuming that only knocked out nucleon 
interacts with the target and that the remaining {A — 1) nucleons of the ini- 
tial nucleus continue to exist in the form of a fragment nucleus (the so-called 
spectator fragment), one can calculate the distributions of these fragments. 



The nucleus consisting of A nucleons with total 4-momentum is described 
by means of the relativistic wave function ^^p^ {[x\'^\ Pi^±]) in which the "longi- 
tudinal motion" of the constituents is parametrized in terms of scale-invariant 
variables 



X. 



(A) 



Po +P3 



Pa,o + Pa,3 

where Pj^^ {ji = 0,1,2,3 is the Lorentz index) and Pa,^ are the individual 
4-momentum of the i-th nucleon in the composite system and the total 4- 
momentum of the composite system as a whole. Square brackets in the argu- 
ment of the wave function denote the set of the corresponding variables 

^"^^ and pi,_L, which satisfy the following conditions: 



X, 



E 



(A) 



< xS^^ < 1; 



A 



Pa 



(5.17) 



The superscript of the variables x^^^^ means that this variable is defined in a 
system of particles whose number is equal to this index. 

The distribution of the spectator fragments in the process of nucleon knock 
out from the nucleus in the laboratory frame (incoming nucleus A moves along 
the z-axis and the target proton is at rest) looks as follows 



E 



da A^/^(sArAr, m^, m^) 



sp 



dP. 



sp 



Ai/2(s,Ml,m2 



a 



NN 



(snn) 



H-^sp, Pi 



sp,. 



aX 



sp 



(5.18) 



Here s is the usual Mandelstam variable for the system consisting of the in- 
cident nucleus A and the target proton, snn is the analogous variable for 
subsystem consisting of the interacting nucleon of nucleus A and the target 
proton. Energy-momentum conservation leads to the following relation be- 
tween these variables 



Snn = s(l - X 



sp, 



p2 



+ 

' sp 



sp 



X 



(5.19) 



sp 



70 



<^NNi^NN) is the total elastic cross section for interaction of the nucleon from 
nucleus A with target. X{x,y,z) is the flux factor, which is defined as 

y,z) = {x-y- zf - Ayz, 

m is the nucleon mass. Ma is the mass of the incident nucleus, Mgp is the 
mass of the spectator fragment, a is given by 

Tfl 

The variable X^p is defined in the following way 

- M + Ea + Pa,s ' ^'-''^ 

Pa,3, Ea and Psp,3; Egp are the z-components of the momenta and the energies 
of the incident nucleus A and of the spectator-fragment (A — 1), respectively. 

The overlap integral /(X^p, -Psp,±) of the relativistic wave functions of the 
incident nucleus and of the spectator fragment is given by: 



.1 ^-1 rfj^-i)' / ^-1 

i{Xsp,Psp,±)= U^A^^^-T.^ 

^ i=l \ i=l 



{A~iy 



x$;(^-^)'([xS^-^)'ri^?' - xr'yp.p,^^bi\^l^,P?])- (5.22) 

The overlap integral is a direct analogue of the corresponding notion which 
appears in the nonrelativistic theory of nuclear reactions (see, e.g., Ref. [39] 
and the references therein). 

The arguments of the wave function of the incident nucleus are related 
to the integration variables and the observable quantities Xgp and Psp,± in the 
following way: 



JA)_^ JA-iy_ Mi)_MiY. „-_io 4_i 

Xi — ay^spXi , P± —p± , I — I, /i i, (tx.o'i^ 

(^) _ 1 Y -(^) _ p ^ ' 

Xa — -L ~ OiA-sp, p j_ — —r^sp,±- 



Thus observation of the spectator fragment makes it possible to obtain infor- 
mation about the character of the "longitudinal" and transverse momentum 
distributions of the nucleons in the incident nucleus A. 

Note that the variable Xgp is defined in an arbitrary Lorentz frame, in which 
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the nucleus A and the target colhde along the 2;-axis as follows 



Y — + Psp,3) /p. ^ .X 

{E^ + P^,3) + {Ea + Pa,3) ' ^ ^ ^ 

where Ps's and E^s are the longitudinal momenta and the energies of the 
corresponding particles. As is easily seen, Xgp is a Lorentz-invariant and scale- 
invariant variable. In the proton rest frame it turns to the form given by Eq. 

mn. 

In the case A = 2, which corresponds to the process of the direct break-up 
of the deuteron, the overlap integral in Eq. flS.lSp is replaced by the deuteron 
relativistic wave function. 



If together with the knock-out of the nucleon from the nucleus the particle 
production takes place, the differential cross section of this process looks like 
Eq. (I5A8|) . but the cross section c^j^ of elastic scattering of the active nucleon 
from incident nucleus on the target nucleon is replaced by the cross section 
cr^^ of particle production in the nucleon-nucleon scattering. 

In Eq. fl5.22p the wave functions <l>j-^'' and "^"/^^^^ of the incoming and outgo- 
ing nuclei are defined in the reference frame, in which their total transverse 
momenta are zero. They are related to the wave functions with arbitrary total 
4-momentum P as follows ^0] 

^^p\[x^\v?]) = 41U([#\P? - (5.25) 



As an example, we consider the following simplest parametrization of the 
relativistic wave functions of the incident and fragment nuclei: 

\ i=i xl I 

and similarly for <I)(^~^) with the substitution A — >■ (A — 1). 

In Eq. f l5.26p is an adjustable parameter, Ca is a normalization factor. If 
the scale-invariant parametrization of "longitudinal motion" of the constituent 
put into the wave function $^'^)([Xj-'^'',p||^'']) is really valid, then the parameters 
a\ would be the same at different energies of incoming nuclei. 

Since we do not distinguish protons from neutrons in our consideration, the 
wave function $'^^^([a;j-'^'',p||^'']) is a symmetric function of its arguments. Solv- 
ing the problem of a conditional extremum under the conditions (I5.17p . we 
find the the wave function f l5.26p has a maximum at zero values of the trans- 
verse momenta of the constituent nucleons and at values of the variables xf"^ 
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equal to 



E mi 

i=l 



Taking into account the relation between the variables x-^-* and Xsp (see Eqs. 
fl5.23p ). we find that the Xgp distributions of the spectator fragments must 
have a maximum at 

Xsp = , ^ . (5.28) 

Note that properties such as the scale invariance of the relativistic wave func- 
tions and the position of the maximum in the Xsp distribution of the spectator 
fragments, do not depend on the concrete form of f l5.26p of the wave functions 
]) and remain valid for their arbitrary parametrization. 

In order to normalize the relativistic wave functions correctly, one has to know 
in general the form of all the interactions inside the relativistic system. Assum- 
ing, however, that the total interaction kernel does not depend on the total 
4-momentum of the composite system, one gets the following normalization 
condition: 

/ n ^ ^ (1 - E -.'^') / n f - E pi') 

•^0 i=l X\ ' \ i=l I i=l \ i=l J 

x|41H[#\p?])r = 2(4vr)^-^ (5.29) 



Substituting the wave function (^^^\[x\^\p^l^]) in the form of (15.261) into the 
normalization condition (15.291) . we obtain the following approximate expres- 
sion for the normalization factor Ca'- 



Ca = y2(47r) V 



■ A 
. i=l 



3(A-1) + 1 

4 



X 



'2af 

TV 



3(A-1) 



exp 



- A 

n 

, i=l 
\ 2-1 



\ i=i 



rrii 



(5.30) 



One of the guiding points in the choice of relativistic wave functions is their 
correct nonrelativistic limit. The nonrelativistic wave function ^'nrUp^), the 
relativistic analogue of which is given by Eq. (15.261) has the Gaussian form 

3(A-1) 

$!^)([p»]) = l^^j exp (^a^EpW^j (5.31) 
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and is normalized by the condition 

A /A 



■' i=i \i=i / 



(A), 
NR\ 



1. 



(5.32) 



Let us expand the combination (pl^"*^ + m|)/a;|"^^ in the power series of the 
parameters p3^/mj, /mi, which are small in the nonrelativistic limit. Re- 
stricting ourselves to the quadratic terms and writing the result in the nucleus 
rest frame Pa = 0, we get 



-(«)2 , 2 



X. 



{A) 



P 



2^ mi + mJ 2^ — ) . (5.33) 



. i=l 



2m 



« \ i=l 



i=l 



mi 



Substituting this expression into the wave function fl5.26p we obtain in the 
nonrelativistic limit (without taking into account the normalization factor) 



^'^\[^\'\p?]) - exp f - a«E "^^^^ ^ 



exp 



2n 



^a[ E"^* 
\ 1=1 



exp 



{A) / 

i=i A J 

^( A \ / A -.{i)2^ 



. i=l 



. i=l 



mi 



(5.34) 



The condition of the correct nonrelativistic limit gives the following relation 
between the parameters of the relativistic and nonrelativistic wave functions 



mi 



A 
i=l 



A 



(5.35) 



For the normalized wave functions, the nonrelativistic limit looks as follows 
^'^r'H^^Kp?]) - v^(2vrW)^Ai$S,^l([j;«]). (5.36) 

Here is the relativistic wave function fl5.26p normalized by the condition 
f l5.29p and ^^^r is the nonrelativistic wave function (I5.3ip normalized by the 
condition (15.320 . Masses of all nucleons are assumes to be equal m^ = m. 

Substituting now the wave functions ^^^^ and <l>("^~^) into the overlap inte- 
gral (15.221) . taking into account the relation (15.231) and integrating over the 
transverse momenta, we obtain 



A-2 



sp,. 



CaCa- 



74 



X exp 



where 



1 - aX. 
X exp 



exp 



Rp2 



aX.Jl - aX. 



sp J 1 



„1 A-1 / A-1 

•^0 i=l \ i=l 



R 

M-1 



X. 



E 



^XspJ U xl^-'\ 



(5.37) 



(5.38) 



Thus, the Psp,_L-distribution is obtained in the analytic form. The integral over 
the variables x^'^'^^ can be calculated approximately by means of the multidi- 
mensional saddle point method [169] in the form of an asymptotic expansion 
in inverse powers of the large parameter a^_i + a^/aXgp. The leading term 
of this expansion has the form: 



MX, 



sp J 



A-l \ 1/4 / A-1 

. i=i / \ i=i 



3(A-2) + l 



TT 



A-2 

2 



a A-1 + a-A/iaXsp)^ 



X exp 



R 



,R 



+ 



aXsp^ 



■A-1 

E 

, 1=1 



(5.39) 



One can estimate the first order correction to the value of integral (15.38^ . 
With account of this correction the integral over the variables x\'^~^^ takes the 
following form 

JiXsp) = UXsp)K{X,p), (5.40) 

where 
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sp J 
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U{A-1) + 



(5.41) 



Substituting into Eq. (15.371) the expressions for the normalization factors Ca 
and Ca-1 from (15.301) and the expression for the integral J{Xsp), we obtain 
finally for the overlap integral 



I{Xsp, Psp,±) = 2(47r 

3(A-1) + 1 



3(A-1) 



A-1 
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(5.42) 



It can be seen from the overlap integral fl5.42p that the distribution with 
respect to the transverse momentum Psp,± of the spectator fragment must 
have a Gaussian form and X^p-distribution will have a maximum at the point 
predicted by Eq. 

Now Eq. (15.181) for the differential cross section of the nucleon knock-out 
from the relativistic nucleus with the overlap integral (15.421) can be used for a 
comparison with experimental data. 



4 Scaling Properties of the Relativistic Wave Functions of ^He and 
^He Nuclei 

To verify the scaling properties of the relativistic wave functions and to extract 
the information about the values of their parameters, we used the experimental 
data on "^Hep interaction at 8,56 GeV/c and 13,5 GeV/c momenta of incident 
'^He nucleus and ^Hep interaction at 13,5 GEV/c momentum of incident ^He 
nucleus. These data were obtained in the 100-cm hydrogen bubble chamber 
|166] , [167j of the JINR at Dubna. We also used the data on the momentum 
distribution of the ^He nuclei emitted at angle 0°65' in the '^Hep — *• ^HeX 
reaction with an incident ^He nucleus having momentum 6,85 GeV/c |168] . 

The experimental da/dXgp and da/dPsp^± distributions of the and ^He 
spectator fragments in the ^Hep — > ^Hpp and ^Hep ^Hepn reactions 
were analyzed in the rest frame of the target proton. (The spectator nuclei 
and ^He were defined as the fragments having the smallest momentum 
among the reaction products in the ^He rest frame.) They are related to the 
invariant differential cross section (15.181) as follows 




(5.43b) 



(5.43a) 



(5.43c) 
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The limits of integration Xsp^-^ and Xgp^^^ were taken from the corresponding 
experimental Xgp-distributions of the spectator fragments. The upper limit 
-Psp,±max is the kinematical limit determined by the condition of the positivity 
of the factor X^sj^n^w? ^rn?) in Eq. flS.lSp : 

= {sX,p - Mlf{l - Xsp) - Am'X^p . (5.44) 

The total elastic cross section cr^^ can be assumed to be effectively constant 
in the considered range of energies and approximately equal to 24 mb |157j . 

The overlap integral (15.421) depends weakly on the values of the parameter 
of the spectator nucleus. Therefore, when the parameters of the relativistic 
wave functions of the final nuclei and ^He are determined by fitting the 
experimental data, there are large errors. For this reason, the fitting was done 
for a fixed value of the parameter of the wave function of the three-nucleon 
nucleus. 

For the determination of the parameter af we used relation (15.351) between 
the parameters of relativistic and nonrelativistic wave functions. It should 
be noted that the parameters of the nonrelativistic Gaussian wave functions 
of '^He, ^He and nuclei are not determined well enough and we do not 
have reliable information about them. The values of the parameters a^^ and 
a^^, given in the literature, vary in rather wide range (see, g.e., Refs. |170j - 
[173] ). The greatest number of data are available for the parameters of the 
nonrelativistic '^He wave function, but the values of the parameter a^^ of the 
Gaussian parametrization (15.311) vary in the range a^^ = 20 — 28(GeV/c)~^. 

The experimental and theoretical Psp,± and Xsp-distributions of the spectator 
nucleus in the ^Hep ^Hpp reaction at 8,56 GeV/c momentum of 
the incident '^He are given in Fig. [191 The same distributions of the ^He 
spectator fragments in the '^Hep ^Hepn reaction at the same momentum 
are presented in Fig. [201 The Psp,± and X^p distributions of fragments in 
^Hep ^Hpp reaction at 13,5GeV/c momentum of incident '^He are given 
in Fig. [21] and momentum spectrum of ^He fragments scattered on 0°65' in the 
laboratory frame in the ^Hep — > ^HeX reaction at 6,85 GeV/c momentum 
of incident '^He nucleus is shown in Fig. [221 The solid curves in Figs. [T9l - 
[22] correspond to the theoretical calculations with values of the parameters 
obtained by fitting the experimental data and presented in Table [H As can be 
seen from this Table the values of parameter af of the relativistic wave function 
obtained by fitting the data at different momenta of the incident nucleus are 
close to each other and agree satisfactorily with the values ~ 5—7 (GeV/ c)^. 
predicted by Eq. (15.351) . 

The positions of the maxima in the Xgp of spectator fragments agree well with 
the values of X^p predicted by Eq. (15.281) . 
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Fig. 19. The Psp,± and Xgp distributions of the spectator fragment in the 
^Hep ^Hpp reaction at '^He momentum 8,56 GeV/c. 




Fig. 20. The Psp,± and Xgp distributions of the spectator fragment ^He in the 
^Hep ^Hepp reaction at ^He momentum 8,56 GeV/c. 




Fig. 21. The Psp,± and Xgp distributions of the spectator fragment in the 
^Hep —> ^Hpp reaction at ^He momentum 13,5 GeV/c. 
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Fig. 22. Momentum Spectrum of fragments emitted at angle 0°65' in the 
^Hep ^HeX reaction at ^He momentum 6,85 GeV/c. 

The approximate energy-independence of the values of the parameters of the 
relativistic wave function in the considered region of energies of the incident 
"^He nucleus allows to assume that in the wave function (15.261) there is no 
other dependence on energy apart from the dependence on the scale-invariant 
variables x^^^ 

For the experimental determination of the parameter af of the relativistic 
wave function of the three-nucleon nucleus, we consider the ^Hep ^ Dpp 
reaction. In this case, the overlap integral has the following form 



X exp 



IiXsp,Psp,±) = 47rcz53^/^m(af)3/2 

^m2(3«X,p-2)2 



aXsp{l 



aX 



sp J 



exp 



^ exp(-4m 



«f + 



aXsp{l 



aX 



spj 



i=l 



aXsp^ 



(5.45) 



Here Cd is a normalization factor 



C 



D 



2^/27^-1/4^1/2, 



5 5 AR AR 



exp[-4m2(af + a| 



-1/2 
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Table 2. The Parameters of Relativistic Wave Functions of Three and 

Four-Nucleon Nuclei 



Reaction, 
momentum 


a^, (GeV/c) 


, (GeV/c) ^ 


2 

X 


/at 


^Hep — > pp 

p.. — Q T^a PpV /p 
-iJJe — ct, (JU vjc V / C 


o - ii-A.eu. 


/ , oy m u, 


da 1 dXsp 

01,10/ lO 


da/dPsp,± 

97 P,'\ /99 


^Hep —>■ ^Hepp 

Pu — S ^fi ripV/r 






54,63/18 


43,15/22 


^Hep — > pp 
PHe = 13, 5 GeV/c 


8 - fixed 


6, 23 ±0,25 


29,01/14 


37,14/17 


Pffe = 6, 85 GeV/c 


8 - fixed 


7,21 ±0,34 


34,93/20 



and and af are the parameters of the relativistic wave function of deuteron: 



i=l 



— a. 



xil — x) 



(5.46) 



which we chose in the form of the relativistic analog of the Garthenhaus- 
Moravcsik wave function 

^r{p) = t^rE^r exp(-afV^). (5.47) 
1=1 



The parameters and are related to the parameters of the nonrelativistic 
wave function (15.471) in the following way 

= A^eMm^ar), af = \ a^. (5.48) 



When comparing the theoretical calculations with the experimental data on 
the distribution of the spectator deuterons in the reaction ^Hep Dpp, 
the parameters of the deuteron wave function (I5.46P were fixed by the val- 
ues obtained in accordance with the expression (15.481) from the values of the 
parameters of the nonrelativistic wave function (I5.47P given in Ref. |174j . 

In Fig. [23] the X^p-distribution of spectator deuterons in the ^Hep — ^ Dpp 
reaction at 13,5 GeV/c momentum of incident ^He is presented. 
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Fig. 23. The Xgp distribution of spectator deuterons on the reaction ^Hep Dpp 
at ^He momentum 13,5 GeV/c. 

The theoretical curve corresponds to the calculation with the overlap integral 
f l5.45p with the value of parameter = 8, 28 ± 0,71 (GeV/c)~^, obtained 
by fitting the experimental data. This value of the parameter agrees well 
with the value af = 6 (GeV/c)^^, predicted by Eq. fl5.35p and used by us in 
previous calculations. 

It is interesting to compare the theoretical calculations on the distributions 
of the parameter spectator deuterons with experimental data at different en- 
ergies of incident ^He nucleus in order to check more precisely the scaling 
properties of the relativistic wave function of ^He nucleus. It is also of interest 
to consider wave functions of a more complicated form. However, our simplest 
parametrizations of the relativistic nuclear wave functions make it possible 
to establish a number of features (such as scale-invariant character of "longi- 
tudinal" motion of nucleons inside the nucleus) that do not depend on their 
specific parametrization. 
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